NPS ARCHIVE 
1967 
GOBBEL, J. 


AN AID FOR THE OPTIMUM DESIGN OF ANTENNA 
SYSTEMS FOR VERY LOW FREQUENCIES 


JAMES THOMAS GOBBEL. JR. 


AN AID FOR THE 
OPTIMUM DESIGN OF ANTENNA SYSTEMS 


FOR VERY LOW FREQUENCIES 


by 


James Thomas Gobbel, Jr. 
Lieutenant, Unit#d States Navy 


B.S.E.E., Duke University, 1962 


Submitted in partial fulfillment of the 
requirements for the degree of 
MASTER OF SCIENCE IN ELECTRICAL ENGINEERING 
from the 

NAVAL POSTGRADUATE SCHOOL 


December 1967 


ABSTRACT 


This thesis will be concerned with the antenna systems used for 
high power radio transmitters operating at Very Low Frequencies (10-30 
kHz). General problems involved in antenna systems, from the viewpoint 
of their equivalent electrical circuits, will be discussed. The re- 
quired matching networks will be covered and their relative performances 
indicated. A brief survey of signalling methods commonly used at VLF 
is included, and the relation between these methods and the aerial 
systems is discussed. 

A computer program (VLFANT) will be developed for a typical 
VLF aerial installation which may be used in the selection of component 
values for optimum operation. Written in FORTRAN IV, it is for use on 
a high-speed digital computer. Circuit parameter values may be changed 
at will, and the system response using different network arrangements 
easily determined. With this program as a design tool proposed systems 
may be completely evaluated and the proposed changes to operational 


systems easily studied. 
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Tx INTRODUCTION 

This thesis will be concerned with the design of antenna systems 
used for radio communications at Very Low Frequencies (10-30 kHz). Prin- 
ciples of design and operation of the systems, from the viewpoint of 
their equivalent electrical circuits, will be discussed. The design 
and use of matching systems will be covered, and a brief analysis of 
information transmission methods used at VLF will be included. 

Selection of the proper complete antenna system is difficult, 
and under field conditions the best designed system may not function 
efficiently. In order to choose component values for optimum opera- 
tion, a computer program (VLFANT) will be developed. Parameter values 
may be varied at will and their effect on overall operation quickly 
determined. 

The computer program is written in FORTRAN IV language for use 
on a high-speed digital computer such as the IBM System/360. The 
U. S. Navy VLF installation at Northwest Cape, Australia has been 
selected as a model and its performance will be discussed. Modifica- 
tion for other installations may be easily accomplished. 

Historically, the Very Low Frequency spectrum was used as the 
primary long range communication medium, with a number of stations being 
built in the period from 1910 to 1912 for world-wide en Begin- 
ning in the 1920's, however, with the discovery and subsequent high 
usage of the "shortwave" frequencies, the utilization of the 10-30 kHz 
range for long distance communications rapidly declined. But within the 
last decade VLF has again come into wide use by both the civilian and 


* Numbers in superscript refer to references, listed at the end of this 
thesis. 
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military communities for accurate time and frequency standards, naviga- 
tion and specialized communications. This renewed interest is due to 
the predictable propagation characteristics of these very low frequen- 
cies, the long ranges attainable, and the ability to penetrate the 
ocean depths, a requirement raised by man's new interest in the under- 
sea pontdes As of July, 1963 over 90 VLF transmitters were known to be 
operating below 30 kHz and in the past few years others have been 
eae The lowest operating frequency regularly in use is 10.2 kHz, 
which is part of the U. S. Navy Omega navigation system, while the 
lowest frequency used for communications is 15.5 kHz at the U. S. Navy 
VLF station at Northwest Cape, Aus Coe 

The earth and lower edge of the ionosphere tend to act as 
concentric shells which guide VLF radio waves. There is little pene- 
tration into the ionosphere by the skywave portion of the signals and 
very small attenuation of the groundwave. Thus signals of notable 
strength may be produced at distances of up to 1,000 km. Absorption 
and attenuation tend to decrease in winter, at night, and with decrease 
in frequency. The major loss in signal strength is a function of 
distance from the transmitter; fading is never gieenrcde Solar acti- 
vity can cause a change in VLF propagation, but the lower the frequency 
the less the effect. For the same reason VLF is relatively immune to 
nuclear iaekout.- ” 

Despite these advantages VLF is not completely utilized because 
the aerial systems are large and inefficient, with much of the generated 
radio frequency power going into ground and copper losses. Additionally, 


the total number of stations which can operate is limited and each sta- 


tion is therefore restricted to a narrow modulation band. 
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A typical modern VLF communications transmitting installation 
uses very high power, often in the hundreds of kilowatts. The antenna 
is very large physically, and the associated components must have good 
mechanical strength and be able to withstand high voltage and current 
guadwenveye In fact, the rf feed-through bushings for Northwest Cape 
are about 16 feet tall and weigh almost 3.5 tons apiece. 

With such an investment in materials, cut and try methods of 
design cannot be used. A proposed system must be painstakingly de- 
tailed, and such variables as stray capacitance included in the 
theoretical analysis. In the past, prior experience with operating 
installations has been combined with modelling studies in an attempt 
to get the best design. Exact mathematical descriptions of operations 
may have been available, but the evaluation of the expressions could 
take months to complete, and with a slight change in a single parameter 
could have been rendered useless. Also the form of a system might not 
fit the nice geometric configurations of a model, calling for approxima- 
tion techniques whose accuracy would be "ball park" at best. 

Today this hurdle is being overcome with the use of high-speed 
digital computers. Commercial antenna manufacturers are designing 
their systems with the aid of computer programs, still retaining models, 
however, as a check on the theoretical results. As one antenna engineer 
has said: 

It is significant that the properties of an antenna 

configuration, for example, can be evaluated more econom- 
ically by computer techniques than by experimental methods; 
however, it is perhaps even more important that the computer 
can also be programmed to change the antenna configuration 
in a systematic way to find a structure that is optimum in 
certain respects. It is not at all unrealistic to expect 


that an optimization that might take several weeks if done 
experimentally can be done within ten minutes by computer. 


ibs: 


This thesis, then, will be concerned with one part of a VLF 
transmitting facility: the interconnection of a previously designed 
VLF antenna with the transmitter. A review of basic principles of such 
a connection will be presented first, and then a computer program 


developed to aid in the optimum selection of the coupling network. 
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Tie BASIC PRINCIPLES 

This section will discuss some of the basic principles appli- 
cable not only to the design of the VLF aerial itself but to the 
matching networks and signalling methods which must be considered in 


the design of the overall antenna system. 


A. Aerials 

It has been stated many times that the most important part of 
a transmitter installation is the antenna system, since the best trans- 
mitter is practically worthless if its energy cannot be properly 
matched to the atmosphere and radiated into space. This is especially 
true in the Very Low Frequency range. 

Wavelengths range from 10 km (16 statute miles) at 30 kHz to 
30 km (48 statute miles) at 10 kHz. The simple half-wave dipole which 
is so popular in the High Frequency bands might be possible here if the 
required real estate was available. But for long range propagation of 
VLF signals, vertical polarization is necessary, and the quarter-wave 
vertical antenna would have to be 7.5 km (12 statute miles) high. 
Consideration of supports and protection against the elements clearly 
makes such an installation impossible. 

Since VLF signals are generally "broadcast" in an omnidirec- 
tional pattern, some form of single vertical radiator seems the most 
appropriate aerial to use. Towers are generally limited to about 
1,000 feet (305 meters), although there are cases where geographical 
features may be utilized to increase the Retene With a limit to the 
physical height of the aerial, it is proper to consider possible 


increases to the electrical height. 
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Consider a transmission line which is one-quarter wavelength 
long electrically, fed with a sinusoidual input, and terminated in a 
resistive load much greater than its characteristic impedance. The 
current through the load will be very small, the voltage across the 
load high. As we move from the load towards the generator the current 
will increase while the voltage will decrease, and a phase difference 
will develop between the two. At the generator end of: the line the 
voltage and current will. be in phase again. If the impedance at any 
point along the line is defined to be the quotient of the voltage 
divided by the current, then we see that due to the difference in phase, 
the impedance will have a magnitude and a phase. With current leading 
the voltage the impedance phase will be negative, the same condition 
as exists for a simple circuit composed of a resistor and capacitor. 

The yertical aerial may be compared to a transmission line, and 
since it is less than a quarter wavelength long, it, too, will exhibit 
a capacitive reactance at the "generator" terminals. If some additional 
capacitance was added at. the top. of the aerial it would appear electri- 
cally to be just more "transmission line" and thus appear longer elec- 
trically. The 1,000 foot tower is so very much shorter than a quarter 
wavelength that in general the more capacitance we can add the longer 
the aerial will appear. This addition of capacitance i known as 
"top loading" and the capacitor itself is known as the "top hat" due 
to its physical arrangement. 

Actual construction often consists of a number of horizontal 
wires run radially from the central tower to outlying towers, for 
distances up to 8,000 feet, forming a circular electrical "plate," 


covering up to two square miles. The other "plate" for the capacitor 
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is the ground, which must be more than just the earth. The basic 
operation of a quarter-wave vertical aerial depends on an 'image" 
aerial of equal length to be formed in the "ground," so the earth 
ground system must form a "ground plane." 

As might be assumed from the discussion above, the proper 
design of a VLF aerial system is relatively complex. Although the use 
of modern analysis has been suggested, it has been the practice up to 
this time to apply basic and empirical formulae and then make a model 
of the proposed erence 

An example of such an approach is the design of the aerial for 
the U. S. Navy station at Cutler, Maine. First, a number of designs 
were studied in theory and many rejected; the most promosing were then 


modeled using a scale of 250071.7+ The best of these models were then 


selected and new models were made with a scale of fev A final 
report was then issued with a recommended antenna, including extensive 
cost data, preliminary sketches, computational results, modelling 
results, and predicted performance against the elements, including high 
winds and ice leadaneee 

Two important parameters for the resulting aerial system are 
the radiation resistance and the effective height, which is "the 
height of the antenna center of radiation above the effective ground 
ver.” The latter is usually determined from modelling measurements. 
The radiation resistance is the true "load" of the aerial and power 
dissipated in this "load" is actually radiated into space. These two 


parameters are proportionally related through a simple expression. The 


radiation resistance normally has a magnitude on the order of one ohm. 


iy, 


Figure 1 shows a sketch of the aerial system in use at the VLF station 
at Northwest Cape, Australia. 

Once the final aerial system design has been made, it is pos- 
sible to consider the system as a series electrical circuit, such as 
shown in Figure 2. . The radiation resistance is the "load," ohmic 
resistances of the tower and cables are lumped into the loss resistance 
term, the capacitance is primarily that of the top-hat, and the down- 
leads to the transmitter form the inductance. 

A more thorough analysis would show other reactive elements 
such as stray capacitance, but they have been combined where possible 
and ignored when practical. It is also possible to use more than one 
downlead but that will not be considered here.” 

In standard network analysis the reactances would be combined 
and the sum considered as a single parameter. This. will be done later 
when it comes to impedance matching, but for the aerial itself the 
reactances must be considered separately. For a typical VLF aerial 
the capacitive reactance can be assumed to be about four times greater 
in magnitude than the downlead reactance. 

The static capacitance (neglecting the downlead inductance) 
determines the major limits for a VLF aerial. The "Q" or quality 
factor is defined to be the ratio between the stored and dissipated 
energy, or the capacitive reactance divided by the radfation resistance. 
The Q algo determines the aerial circuit bandwidth. For VLF aerials 
the Q has a value of several hundred, and the bandwidth is measured in 


ten's of Hertz. 
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where Xx. = capacitive reactance at the 


operating frequency 


ie = operating frequency 


There is a second set of values for the aerial Q and bandwidth 
often encountered which involves the loss resistance. At the input 
terminals of the circuit in Figure 2 the loss and radiation resistances 
are seen as a single resistance sum. From the network analysis view- 


point, Equation (1) must be modified as follows: 


jx, 


Q a ngs 


one Readeetien a5 ices 

If an aerial system with the loss and radiation resistances 
equal is considered, the Q as seen from the input terminals is only 
one half that of the aerial itself and the bandwidth is doubled. Since 
the losses depend upon the actual aerial installation and the associated 
coils, leads, ground, etc., the Q in Equation (1) is the one of most 
concern in the initial aerial design. From the operational viewpoint, 
however, the Q determined from Equation (1l.a) and the corresponding 
bandwidth are the parameters of most concern. 

One of the two major VLF system limitations appears here. This 
is the breakdown voltage rating of the various components in the aerial. 
(The second is a limitation on the power the transmitter can produce.) 
With the low radiation resistance, high power necessitates high current 


which in turn produces high voltage across the reactances. The voltages 
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Figure 2. Basic Aerial Circuit 
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Figure 3. Voltages and Power 
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across the capacitive reactance (Figure 3) and the total aerial react- 
ance are the “top-hat" and "base" voltages respectively. Note that the 
voltages may be expressed directly in terms of the quality factor for 


the respective reactances, when the aerial is tuned to resonance. 


Segoe tee = One ae ons) 
\yene F Sdeantlesd Vn ~.- (4) 
where V5 = voltage applied to the 
a series circuit 
Qaownlead 7 2 Of the downlead or total reactance, 


defined similarly to Equation (1) 

Another problem is corona which takes place just before break- 
down occurs. Corona develops at high voltage points which are physi- 
cally near points at ground potential, and it may occur at the endpoints 
of the top-hat, the base of the downlead, cable connections, sharp bends, 
and at points where the wires may be frayed or nicked. High rf noise is 
generated and the aerial circuit is detuned. Corona rings and the caging 
or bundling of cables may be used for protection. 

The aerial self-resonant frequency is the series resonant con- 
dition when the inductive and capacitive reactances are equal and is 
usually outside the VLF band. The top-hat provides the major part of 
the aerial capacitance, without it the self-resonant frequency would be 


even higher. 


1 
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As the aerial capacitance is increased the Q will decrease and 
the bandwidth increase. The effective height will also increase and 


with it the radiation Seegigeraa” The top-hat voltage (a function of 


Ze 


radiation resistance and capacitive reactance) will tend to decrease. 
Thus the greater the top-hat capacitance the better the aerial will 
perform. 

The aerial efficiency is the ratio of the radiation resistance 
to the sum of all resistances in the circuit. Since the radiation 
resistance is of the order of one ohm, all loss resistances must be 
equal to or less than this to obtain an efficiency of 50% or greater. 
In practice this can be done by using large diameter cables to reduce 


losses. 


Arad iat ion 
Efficiency = <i 0) 


é + 
radiation Eioee 


Equation (6) is good only at a specific frequency since the radiation 
resistance will normally increase with frequency. The lower the oper- 
ating frequency, the more carefully all losses must be considered. 

In summary, the normal VLF aerial is much shorter than an 
electrical quarter wavelength and is electrically lengthened through 
the use of top-loading capacitance. It is physically large (mainly 
due to the top-hat), has a high Q and narrow bandwidth. The achieved 
efficiency will in practice be close to 50%, and limitations on power 
handling ability is usually caused by breakdown voltages across the 


reactive elements. 
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Be Impedance Matching Networks 


Maximum power transfer between a generator and its load occurs 
when the load resistance is equal to the generator resistance. If the 
generator impedance has a reactive component then the load should have 
a reactance of equal value but opposite sign. The generator and load 
impedances are then "conjugates." 

In radio frequency applications the generator can usually be 
made purely resistive, and thus the load should be a pure resistance. 
But even with a pure resistance load maximum power transfer may not 
occur. For example, if a vertical quarter wavelength radiator, half- 
wave dipole or folded half-wave dipole are each resonant at the same 
frequency their input impedances are purely resistive. But their 
magnitudes are 36, 72, and 300 ohms respectively. If the generator 
has a resistance of 50.ohms, maximum. power transfer cannot be directly 
accomplished with any of these systems. 

A network may be added between the generator and antenna which 
will allow maximum power transfer. The generator sees the network and 
antenna combination as a 50 ohm resistive load, and the antenna looking 
into the network-generator combination sees a pure resistance source 
of the proper magnitude. The network which meets these requirements 
is known as an impedance matching network. 

It is desirable to introduce as little insertion loss as pos- 
sible, so the network should be "lossless.'' If resistors are used in 
the network such a requirement cannot be met, but by using high quality 
reactive components the lossless case can be closely approximated. No 
component is ideal, but the ohmic resistance can be made very small, and 


in any case, much, much less than the alternating current reactance. 
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VLF aerials have been shown to have a resistive component in 
their input impedance, and since inductive and capacitive parameters 
are present, also a reactance term; except at the precise resonant 
frequency. When operation at other than resonance is desired the 
reactance term will appear and increase in magnitude the further 
operation is attempted from resonance. Below resonance the reactance 
will be capacitive; above resonance inductive. In the following discus- 
sion it will therefore be assumed that the aerial is a complex 
impedance load. 

Impedance matching networks are built using the basic prop- 
erties of series and parallel circuits. In Figure 4 two circuits 
with identical behavior at the same frequency are shown, the second 


being the parallel form of the first. The parameters are defined as 


follows: 

Z = R. te 5X, SG) 

as = Re =F i sles) 
where 

z = impedance of series circuit 

a = impedance of parallel circuit 

R, = series resistance term 

jx, = series reactance term 

a = parallel resistance term 

i = parallel reactance term 


The impedance of the two circuits is the same at the same fre- 
quency. The quality factor of "Q" defined in Equation (l.a) for the 
aerial may now be generalized to apply to both the series and parallel 


circuits. At the operating frequency, both circuits have identical 


mors." 


2) 


Q, aor for the series circuit eee) 
Ss 
R 
Q, = TFs] for the parallel circuit veel) 
P 


Equations for converting between series and parallel circuits 


are given below: 


Parameter "'Q'' Equations Reactance Equations 
2 3 Z 
x = R6(Q. + 1) or Ro[1 + Re ] pe ciil) 
Qe ar dl R, D 
ie _ ao or xX [1 + ~ ] euiGl2) 
R R 
R 2 R es) 
= or Res. 
2 
: ee [1 +27) 
P 
jx, 4X 
4X 2 R (14) 
= or 
in D one 
pe [2 +2)*1 
o; P 


Either the "Q" equations or the reactance equations may be used 
in transforming from one circuit to another. When the reactance of a 
series circuit is greater than the resistance by a factor of 10 or meee 
(or less by a factor of 10 for the parallel circuit) then the equations 


above may be simplified. 


2, oe 2 
a = (Q.) R. or re Pa LES) 
ix, x 5X. or 5X. RS) 
R R 
R ie 2, i 

' ~ (Q) or RS 9 olay 

G ) 

P 
5X, ~ ios or 1p mee Gh) 
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Examination of equations (15) through (18) shows that an error 
of 1% occurs when Q is 10; the error decreases as Q increases. 

If the circuits in Figure 4 are to be modified to be purely 
resistive an additional reactance is needed. Thus if a reactance of 
the same magnitude as Xo but of opposite sign is added to the series 
circuit, the reactance will sum to zero and a pure resistance results. 
By the same token, adding a reactance which is the conjugate of - to 
the parallel circuit will give a resultant impedance equal to Be 

The above operations may be combined to produce different 
results. If a reactance of magnitude Xe but of opposite sign than i 
is added in series to the parallel circuit, then the total impedance 
becomes R.- This may be seen by first transforming the parallel cir- 
cuit to its series equivalent and then adding the new reactance. 
Alternately, if a reactance of magnitude with sign opposite Xe is 
added in parallel to the series circuit, RS is the resultant impedance. 


These results are shown in Figure 5 and summarized below: 


Circuit Added Reactance Connect In Resultant R 
Series = )3X, | series R, 
gots z |3%, | parallel R, 
Parallel - )3x, | series R. 
Parallel = |3%,| parallel Bes 


The procedures above are the basis for a common type of imped- 


ance matching networks. Methods for connecting a reactance to a cir- 


cuit and obtaining either the series or parallel value of the resistance 


have been shown. But suppose a resistance value other than RS or 3 is 
desired, but still without any reactive terms present. This can be 


done but may require a second reactance to be used. 


2h 


Figure 4. Basic Series and Parallel Circuits 


mre J. Gastie Matching Circuits 
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Following the work of Fowler and Kemp, the circuit shown in 
Figure 6 is iecdnime With the circuit broken at A looking to the right, 
a resistance Ss is seen. Similarly, if the circuit is broken at C, 
looking to the left R. is seen. Thus resistances R, and _ are matched 
to each other. However, at cut B the circuit to the right is R - JX 
to the left is Sy in parallel with so If the series circuit to the 
right is transformed to its parallel equivalent, Be is in parallel with 
=," The left and right hand circuits are seen to have equal resistive 
components oS but conjugate reactances (equal magnitude but opposite 
signs). The parallel circuit could have been transformed to a series 
circuit and given a resistance of R. at cut B. No matter what form is 
used, the impedances at B are seen to be conjugates. 

From Equation (11) Ss is seen to be always. greater than Ro: 
Then to match two resistances the higher resistance may be connected 
across the shunt arm of the network to the right of cut A. The two 
reactances are adjusted to give conjugate impedances at cut B. If the 
two loads to be matched have reactive as well as resistive terms, then 
the reactances in Figure 6 are first adjusted to make the loads purely 
resistive and then modified to meet the matching requirements. Due to 
its physical shape, the network between cuts A and C in Figure 6 is 
known as an "L" network. 

Figure 7 shows two other forms of matching networks, the Tee 
and the Pi. These networks are seen to be the combination of two L 
sections combined in different ways. They may in fact be designed as 
two L sections and then combined. 

When pure resistances are matched the two reactances in the L 


network must have opposite signs, in the Tee and Pi this is extended 
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“izure 7. Tee and Pi Networks 
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to having one element different from the other two. These restric- 
tions do not apply if the loads have reactive components. 

If the Tee and Pi networks are constructed as shown in Figure 7 
then they are called symmetrical.’ The inductors and capacitors may be 
interchanged in either network without destroying the symmetry. It 
should also be clear from Figure 7 that the Tee and Pi networks shown 
are duals, and at a given frequency one may be directly substituted 
for the other. 

Tables I, II, and III give the equations needed for designing 
impedance matching networks of the L, Tee, and Pi forms respectively. 
In all cases the two impedances to be matched are assumed to be purely 
resistive. If they are not, the networks should be designed on the 
assumption that they are, and then the network legs may be modified 
as required. 

An alternate method of computing values for the Tee and Pi 
networks is to consider each as a lumped parameter equivalent of a 
quarter-wave transmission line. Such a line is often used for imped- 
ance matching, and has the additional property of inverting the load 
as seen at the input network terminals. That is, if a load which has 
infinite impedance (open circuit) is connected to the Tee, reactances 
Xo and XX, are seen to be a series resonant circuit, which has zero 
impedance. Conversely, a zero impedance load (short circuit) causes 
x and XX, to be a parallel-resonant circuit which has an infinite 


impedance. The Pi network may be similarly described. 


oi 


TABLE 1 


L SECTION DESIGN EQUATIONS 


Diagram 


j Assumptions 


Equations 


Keand x are taken with 
a 


opposite signs. 
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TABLE IT 


TEE NETWORK DESIGN EQUATIONS 


Diagram 


€ 
i > 
Assumptions R. R, 
Both R, and R. are purely resistive. 
> 
Ro both R, and R, 


__ 


Pe 721)) 
ae G2) 
Equations 74-9) 
R 
- ro 
Xa = ‘E es) 
Re Z 
2 
R 
Se 
kg se -1 «se (2op 


Note: Xx, and X, have opposite signs. 


d 
Xx and xX, have opposite signs. 
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Given R. ; R, ; Xx 


a 2 
Rk =@eepl +  ] wnG26)) 
c a R 
a 
Re Z 
Xi, 7 X, [1+ ] suns G27DE 
a 
Equations 
x Use equation (22) 
Xx, Use equation (23) 
xX, Use equation (25) 


Xx Use 
a 

} Equations xX Wee 
@ 

Xa Use 

Xx Use 
e 


Equations 


R, 


equation 
equation 
equation 


equation 


opel 
Ra = R, [1+ @ 


] » ++ (28) 


(21) 
(23) 
(24) 


(25) 


Use equation (24) 


Use equation (25) 


Use equation (26) or (28) 


TABLE II (continued) 
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TABLE ITI 


PI NETWORK DESIGN EQUATIONS 


Diagram 


Ry ae 


Assumptions [ Both Ro and RX. are purely resistive. 


< 
| | R. both R. and Ry 


iver | 


| Equations 


Note: Xx. and Xa have opposite signs. 


xX and xX. have opposite signs. 


2) 


Equations 


Equations 


- 


TABLE TL (continued) 


a2 
i ) 
a 


Use equation (32) 


Use equation (33) 


Use equation (35) 


CC MNES es 


1+ ey? 


Use equation (31) 
Use equation (33) 


Use equation (34) 


“RX, ir, + Ae - sie +6 (40) 
“RX, ifr, 2/RR, - me oa (49D 


Use equation (34) 
Use equation (35) 


Use equation (37) or (39) 
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The networks described, although designed for impedance 
matching, will also act as filters. In work at radio frequencies one 
is usually interested in surpressing higher frequency harmonics, and 
thus a low pass filter response is desired. Or there may be undesirable 
sidebands produced by modulation of the basic carrier system, and then 
a bandpass filter which would discriminate against unwanted frequencies 
would be called for. 

The choice of which impedance matching network to use is based 
on a number of factors. What is the ratio of the two resistances to 
be matched? What type of reactance do the loads have, is it large or 
small, and how may it best be cancelled out? Will one of the loads 
vary significantly with frequency or will the resistance remain fairly 
constant and only the reactance vary? How many components are avail- 
able for the network? What power and breakdown voltage ratings must 
be met? What magnitude of stray capacitance will be encountered? 

(The Pi network with capacitors for the legs would be useful here.) 
What about ease of operation; can tuning be quickly. and accurately 
accomplished by the personnel involved? 

Additionally, what power will be dissipated by the network? 

The proper choice of network may mean a substantial improvement in 
efficiency. Terman points out that "for a given [impedance] transforma- 
tion ratio, the "L" section has a lower loss than either the Tee or 


pi"? 


Example: 
As a particular example of impedance matching, the aerial 


circuit shown in Figure 2 and redrawn in Figure 8(A) will be used. The 


impedance values shown are typical of a VLF aerial, and are chosen to 


a 


show the range of magnitudes encountered. It is desirable to use a 
matching network to obtain a purely resistive output of 12.5 ohms at 
the operating frequency. 

A number of variations of this circuit are possible. In 
Figure 8(B) the reactances have been combined, since this is what a 
radio frequency bridge connected to the input terminals would see. 

It is also possible to take the capacitor and resistor as a sub-network 
and convert them to their parallel equivalent circuit, as shown in part 
C. Since the Q of this circuit is greater than 10 by Equation (9), the 
transformation may be made using Equations (11) and (12). A similar 
procedure yields the circuit shown in part D. . The parallel circuit of 
part E is determined from the circuit of part B. 

Since the desired resistance at resonance is 12.5 ohms, and such 
a value does not appear from any. of. the manipulations in Figure 8, the 
most direct approach to the desired result would be to use the circuit 
in part B and add a matching network. Figure 9(A) shows such an approach. 
First a reactance of equal magnitude. but opposite sign is added in series 
to produce a pure resistance. Then a simple L network is designed to 
convert the 0.308 ohm load to the desired 12.5 ohms. When constructed, 
the series inductance would be included in Xx. so that the completed 
network would be to the left of the dashed line. Values for the L 
network are found from Equations (19) and (20). 

A.second approach would be to use the parallel form of the 
aerial circuit. Here the new inductive reactance is connected in 
parallel and the L networks used to convert 7,550 ohms down to 12.5 
ohms. Again the completed L. network would be the circuit to the left 


of the dashed line. 
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Figure 8. Variations of Basic Aerial Circuit 
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The networks in Figure 8, parts C and D may also be used. If 
it is desired to have a pure resistance at resonance the total series 
reactance must have the same magnitude but opposite sign of the parallel 
reactance. In both of these cases it is seen that an additional induc- 
tive reactance of j48 ohms is required and the total network resistance 
will be the same as the series resistance, 0.308 ohms. 

There is also the possibility of using some form of transformer 
arrangement to achieve the desired 12.5 ohms. In fact a single induc- 
tor could be used as an autotransformer if the taps were properly 
chosen. If the circuit in Figure 8(E) has an inductive reactance of 
348 ohms connected in parallel across it the result will be a parallel 
resonant tank circuit with a resistance of 7,550 ohms. If the inductor 
is tapped as shown in Figure 9(C) then 12.5 ohms will be measured 
between the tap and ground. 


The relation. between turns and impedance for an autotransformer 


is given by: 
Z 
vk, = Gy 
i 
where Z, = primary impedance 
Zo = secondary impedance 


N = turns ratio 
With Z4 and Zo equal to 12.5 and 7,550 ohms respectively, N is 
found to be 1/24. A total inductive reactance of 448 ohms is needed, 
and since inductance varies directly with the number of turns, a selec- 
tion of the tap such as that shown in Figure 10 can be made. Here it 
should be noted that the same results may be obtained by considering 


the circuit shown in Figure 8(B) to be made into a series resonant 


Gime. 
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Figure 9. Various Matching Arrangements 
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546 


Figure 10. Autotransformer Matching 


Figure 11,7 simplified Circuit 
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Figure 11 shows the resultant circuit. It is more convenient 


to consider the series form of the aerial circuit, which has an input 


impedance defined Bees 


2 a 
RG) + j([X,R ac one) 


Z = 2 2 
in R + (x, - X.) 


«+ (43) 


At resonance the reactance terms will partially cancel, and 


the simplified form below results: 


(x, ) 
ee eo JX, ones) 


Thus by tapping the coil to lower the input resistance a 


Z 


reactance term has been added. (It should be noted that the proper 
tap on the coil is most easily. found by Equation (44), since the 
resistive part of Zon and the value of R are shown. Xr may be 
computed, and the total coil is merely equal in magnitude to the 
capacitive reactance.) 

Since the reactance is small it seems apparent that at some 
frequency other than.resonance it will become zero. Figure 12 shows 
the general circuit, where the series capacitive reactance is assumed 
greater than the series inductive reactance at resonance. If the 


series sum is denoted by es Sr the input admittance for the circuit 


is given by: 


R i ee _ u 
Yon = R2 +X 2 Rr? +X 2 x aco 
sum sum Lo 


The real part of the admittance at resonance should be 1/12.5 


mhos» 80 solving for ae 


x = fR (12,5 - R) +3 (46) 


suin 
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Figure 12. Zero Reactance Circuit 
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Figure 13. Final Circuit 
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With an R of 0.308 ohms the positive root fur the magnitude 
of amas will be used, which has a value of 1.93 ohms. It is now neces- 
sary that the reactive part of the admittance be zero so solving for 
the parallel inductive reactance: 


R- ees ; 


x - aaa se Gad) 
ie sum 

This inductive reactance is found to be 1.99 ohms, and the 
complete circuit may now be finished as shown in Figure 13. Comparing 
Figure 13 with Figure 10(B), the value of the inductance has been 
increased very slightly, but the total input impedance of the circuit 
at resonance is now purely resistive. 

Examination of the circuit in Figure 13 shows the part to the 
left of the dashed line to be a L network. When values for the L 
network shown in Figure 9(A) are compared with those in the final cir- 
cuit, they are seen to be the same. Thus after some involved mathema- 
tics the same results as those obtained from the comparatively simple 
Equations (19) and (20) are determined.. Several important conclusions 
can be drawn from this, when deaigning a matching network. 

First compare the magnitude of the desired resistance (12.5 
ohms) with the values obtained from manipulating the aerial network. 
As a series circuit the aerial had 0.308 ohms resistance, with the RC 
parallel configuration the resistance was 12,320 ohms, with the RL 
arrangement 715 ohms, and with the reactances summed and then formed 
into a parallel circuit 7,550 ohms. Since the series value of 0.308 
ohms is the closest to 12.5 ohms, use the series circuit. 


Next compute the values for the L, Tee and Pi networks which 


will make the desired transformation. At this time neglect the aerial 
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reactance completely. The L network will probably be the easiest to 
use although the wider tuning range available or size of components may 
show one of the others to be preferable. 

Finally include the aerial reactance (here ~-448 ohms) in the 
circuit and modify the impedance matching. network values to achieve 
the desired values. Choose variable tuning elements for the network 
so that any variation in component values in the final circuit may be 
compensated for. 

One additional word about the circuit shown in Figure 13; it is 
simple and requires only two components, the values are within easily 
obtainable ranges, and the aerial is provided with a direct current 
path to ground which provides good protection against lightning and 


other electrical discharge. 
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C. Transformer Coupling 


A coupling method not mentioned previously involves the use of 
transformers. They may be constructed to meet many different impedance 
ratios, but generally the magnitude of the loads must fall into a 
relatively narrow range. For example, a transformer with an impedance 
ratio of 10:1 and designed to match 500 to 50 ohms probably would not 
work well if used to match 50,000 to 5,000 ohms even though the ratio 
is the ite 

Transformers utilize iron cores at low frequencies, but at radio 
frequencies losses in the core become excessive and air cores are used. 
Satisfactory operation is usually limited to fairly small deviations 
from the center design frequency. 

At radio frequencies transformers find use where a relatively 
constant impedance ratio is involved, only small deviations in fre- 
quency take place, and as little attenuation as possible is desired. 

An additional advantage is that there is no direct electrical con- 
nection between the primary and secondary circuits. Energy transfer 
is made through magnetic coupling, a fact of some importance where 
high direct current voltages such as those in power amplifiers are 
present. 

A transformer wound with low resistance wire may be considered 
as a lossless two-port reciprocal network. If Z,, is defined to be the 


iit 


input impedance, Z the output impedance looking back into the network, 


Ze 


and Z,,, equal to Z,, to be the transfer impedances, then the impedance 


12 Pa 
matrix for the transformer may be written: 
Z Z 
2 ue > + 6 (48) 
ae @2) 
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Considering only the transformer in Figure 14: 


BE, = 2,1, + 42%, » ++ (49) 

E ~ 254 I, ap Zo9 I, wmco0) 
If the transformer is lossless and terminated in an impedance Laesa> 
then: 

E, = - I, 24 oad ys cenlly) 


By manipulation of equations (49), (50), and (51), the input impedance 


looking into the transformer will be: 


ih 
Bee -T. 
aL 22 
(a) 
or ZL = Zz ae ae eS Se e°8 A52) 
in eal Zo9 + ee | 


For the lossless case, assuming no resistance in the input, 
output, or transfer. impedances. (or else the reactance is much greater 


than the resistance), Z, _.may. be written as: 


in 
2 
. - (wM) 
7 a ie os eo 
dn ul jwL, a Z4 oad 
where Ly = primary transformer inductance 
Ly = secondary transformer inductance 


M = mutual inductance 
w = frequency in radians 
A common application of transformer coupling is in the matching 

of a vacuum tube to a load. The tube may be considered as a generator 
in series with a high resistance, while the load is usually a much 
lower impedance. Such a circuit is shown in Figure 15, Examination 
of Equation (53) shows that the transformer-load combination may be 
replaced by a complex impedance (Figure 16). A capacitor is often 


added to the primary circuit to tune out the reactive part of the 


new load, and the desired impedance matching can be accomplished. 


48 


Figure 14. Basic Air-Core Transformer 
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Figure 15. Vacuum Tube Source 
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Figure 16. Equivalent Circuit For Transformer 
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If the vacuum tubes are to be used as power amplifiers, the 
power can be increased to the limits of the tube without problen, 
providing the transformer windings will handle the high voltage and 
current present. Such units can be built, although they become 
physically large in order to keep the losses low and prevent any 
voltage or current breakdowns. 

When one amplifier has reached its power limit additional 
units may be connected. Figure 17 shows two such units. With several 
units on the line, the transformer secondaries and the load form a series 
circuit and the same circulating current will flow through all the series 
elements. Since all transformers are assumed. identical, each will have 
the same voltage rise across the secondary terminals, the sum being 
equal to the voltage drop across the load. Then with four units on 
the line, Figure 18, each will supply one-fourth the voltage across 
the load, and with the same current as the load, in effect see one- 
fourth the impedance of the load, by Ohm's Law. If the transformer 
secondary reactances are cancelled out. at the operating frequency 
by a series capacitance in the circuit, then in effect with multiple 
units each will see a load resistance equal to the original load 
divided by the number of units in use. 


Modifying Equation (53): 


2 
_— (wM) 
ta dwt jw, + Load, + (o@ 
=. cwe + Zload 
Load | = Gay 
where C = secondary series capacitance 
where Load = portion of original load now seen 


by each transformer 


n = number of units on the line 


30 


PUT, 


SSTISS 


i}? SASEIT Tay Teveneone 


"2 Sdusty 


CH JSTITTAuy sero 


| 


Ee 


CLO) 


(a5 
) 


ea maar © —| =v NO ae ae 


DL 


UCTIVdsedGg Selaes 


4 TU; 


AVC w "OT edits Tn 


} O 
CH 
4TUy 
{ 
j 
\ 
WL --- 5 G 
1 
<. => ! 
T# 2TYUn ee ae 
| q 
Ig eee: cay 
ne cS 


ae 


With the secondary load of each transformer now a pure 
resistance, the primary circuit may also be made purely resistive by 
the tuning capacitor. When operation is at other than the normal 
operating frequency the secondary load will become reactive, but so 
will the primary circuit, and since all units are assumed identical 
all will behave in the same manner. 

Two problems now appear. First, and really most important, 
the power amplifier tubes are operated Class C. This means that 
they are on for less than 180 degrees in each cycle, and the plate 
voltage is not a constant. This entire area of Class C operation is 
a very important one which should be thoroughly investigated. Such 
complications can be temporarily avoided here by taking the time- 
average or dc value of the plate voltage and assuming it constant. 

Secondly, there will be some interaction among the power ampli- 
fier units, which are normally tuned up in sequence. As each additional 
unit is brought on the line, tuning conditions for those already present 
will change. Since the units have no magnetic coupling, and are fed 
from the same signal source, the interaction should not be too great 
in practice. However, the description above neglects such problems, 
and it is an area which requires a more detailed analysis. 

The approach above is easy to use: just divide the normal 
load by the number of units on the line to find the load for any unit. 
But the restriction that each unit be identical must be completely met. 
The only true test is to actually compare the units which will be used 
under the operating conditions, for after all, this is the information 


which must be used in actual operation. 
The effects of transformers in matching impedances has been 


discussed above. The term "M," shown in Figure 14 and used in 


a0 


Equation (53) also has an important effect on the energy transferred 
from the primary to the secondary. The emphasis in the following pages 
will be on how this energy is transferred. 

A transformer consisting. of two windings may be used to couple 
energy from one circuit to another, without any direct electrical 
connection between the two. The coupling takes place when an elec- 
tromagnetic flux is common to both, and a mutual inductance is created. 
The flux, number of turns.in the primary and secondary windings, and 
primary and secondary current may be manipulated to determine the 
mutual inductance in henrys. 

A form which is usually of more convenience when working with 


transformers already designed is given below. 


M = k Ly L, Are Se) 
where M = mutual inductance 

Ly = primary inductance 

Ly = secondary inductance 

k = coupling coefficient 


The coupling coefficient is a function of the geometric 
relation between the primary and secondary coils, theoretically equal 
to zero when the coils are infinitely far apart and equal to one when 
the coils are infinitely close together and the flux generated linking 
both coils fully. The mutual inductance M may range from zero to a 
maximum given by Equation (56) when k is one. In actual practice 
k cannot equal one, but may be made very close to one. (The auto- 
transformer discussed in the Impedance Matching section has a coupling 


coefficient of one.) 
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Transformers are frequency-dependent devices, and in most 
applications both the primary and secondary are tuned to the same 
frequency. Under these conditions the effect of coupling variation 
may be seen directly as shown in Figure 19. Four cases are shown with 
the variation of primary current and secondary voltage versus frequency 
sketched for each case. It is also possible to show the relation be- 
tween primary and secondary voltages or currents; the choice depends on 
the particular interests of the designer. 

"Undercoupling"” is defined to occur when both the primary and 
secondary sides of the transformer are tuned to the same frequency and 
an increase in coupling will cause an increase in secondary voltage. 

If the coupling has been chosen such that the maximum secondary 
voltage is realized, then the circuits are said to be "critically 
coupled." 

As coupling is increased beyond that required for critical 
coupling, the secondary response curve is seen to "widen," that is, 
have a small variation in magnitude as the frequency is varied about 
the center frequency. When the widest response is- obtained without 
any peaks on the curve, then the circuit is "transitionally" coupled. 

Coupling may be increased still further and two peaks appear. 
This condition is the "overcoupled" case. 

Transformers operating at radio frequencies have low loss 
windings such that the coils have an extremely high Q and any ef- 
fective resistance in the circuit comes from external parameters. 

Thus the system in Figure 20 will have R, and R, as the dominant 


i 2 


resistances. The capacitors Cc, and Cy are used to tune out the 


inductive reactance of the primary and secondary respectively. 
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Transformer Circuit 


Double-peaked Response 


The Q of the primary circuit may be defined as: 


Q = [Sa = a eee (oy 
il Ry Ry 
and for the secondary circuit: 
Q = —y eee (58) 
Ze R, R, 


where the reactances are measured at the resonant frequency. The 


coupling coefficient for critical coupling may now be defined as: 


1 
k = ... (aap 
oo Aree 
By substituting the expressions for QQ and Q from Equations 


($7) and (58). ko may also be defined as: 


Ri 


2S oss (OOD 
c xX 
ei as 


It was shown earlier that the square of the mutual inductance 
is a factor in reflecting impedance from one side of the transformers 
to the other, so it is of interest to determine the value of Ga 
when there is critical coupling. Thus in Equation (56), using 


critical coupling as defined in Equation (60): 


M = oo COM 
Solving for the mutual inductance reactance squared: 
Z 
(wM)” = Ri R, ose (Oey 
Examination of Equation (62) indicates that for the undercoupled case: 
2 
< 
(wM) R, R, 2s (ap 
while for the transitional and overcoupled cases: 
2 
(wM)~ > R R, oo (GA) 
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In radio-frequency circuits transformers are used for impedance 
matching and ease in coupling power from the transmitter or other source 
to a load. It is desirable to do this as efficiently as possible, so 
there should be maximum secondary power. 

Bandwidth considerations are also important, since sideband 
intelligence should be passed while undesired harmonics are simultane- 
ously attenuated. A compromise between power and bandwidth may be 
required for best overall operation. 

Critical coupling occurs when the secondary voltage or current 
reaches its maximum possible value. At this point the primary response 
has two peaks, the secondary still has one peak. If the Q's for the 
primary and secondary are equal then further coupling will cause the 
primary peaks to become more prominent and peaks will appear in the 
secondary. The transformer is then overcoupled. 

At high frequencies where modulation sidebands are only a small 
percentage of the carrier frequency, the critically coupled case may 
be optimum. At VLF, however, the aerial circuit is very narrow in 
bandwidth and it is desirable to have other circuits with bandwidths 
greater than that for the aerial. Thus it should be determined that 
the transformer will not place any additional restrictions on the 
system response. 

If the primary and secondary Q's are made unequal, then an 
additional coupling response is obtained between the critical and 
overcoupled conditions. The primary response will show peaks as 
critical coupling is reached and exceeded, but the secondary response 


will widen somewhat before the single peak changes into double peaks. 


By) 


It is desirable to determine when the secondary response will 


aly 


change from a single to double peaks. This can be shown to occur when: 
5 ee 
(wM)” = ———— »sa(65) 


If the mutual inductance reactance squared is less than the term on the 
right then there will be only a single peak (undercoupled condition); if 
greater there will be two peaks (overcoupled case). If R, is equal to 


R,, the peak or peaks will be equal to the maximum response at critical 


2? 
coupling [Equations (62) and (65) are then identical]; when the two 
resistances are not equal the peaks for overcoupling will be less than 
for the critical case. 

It is of interest to examine the magnitude for the various 
coupling conditions above. Figure 21 shows a double peaked case, and 
it will be assumed that Ry and Ry are not equal, and the primary and 
secondary Q's are likewise not equal. Following the work of Henney, 


the gain at the center frequency for any coupling coefficient whether 


with single or double peaks may be written: 


Gain = ———__"»—_____ - __™"__ suaselOooy 
wC, C, [R,R, + (wM)~] oancH +k 


This magnitude is marked in Figure 21 by the dimension "c." 


The maximum possible gain at the center frequency occurs with 


critical coupling for either single or double peaked curves. It is 


proportional to dimension "p'' in Figure 21 and is given by: 


X 
al is Oi oN ie 
Gain = ———— a= = oo s (Ome 


2 Ry Ry 
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When the circuit is overcoupled the dimension "p" represents 
gain at the peaks of the two-peaked curve and is given by: 


Gain = $—_—_— ———E———————————— os» GOOD 


| 2 Z Z 2 
we, C, (2wM) (R, + R,) - (Ry - R, ) 


The location for the peaks for the double-peaked circuit 
may be located by using the following equation, where - is the 


resonant frequency, and f the peak escucaey —! 


~» (69) 


The dual signs in the denominator are for the upper and lower fre- 
quencies f1 and f,- The bandwidth between peaks (Ap) may be computed 
from Equation (69), and the bandwidth between points where the gain 
equals that at resonance is given by Ac: 

we <(s Ap 56 (7/0) 

If a transformer has already been constructed there is an 
additional method of determining its coupling coefficient in the field. 
The primary circuit is tuned to resonance at frequency fo with the 
secondary winding open. Then the secondary is shorted and the primary 


re-resonated, this time to the higher frequency fe: The coefficient 


of coupling is then given by: 


ic fF 
k= -& Pann) 
s 


where fy = primary resonant frequency, secondary open 


f = primary resonant frequency, secondary shorted. 
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In summary it may be stated that the peak gain of the secondary 
response curve depends on the ratio of (wo) 2 to the product of Ri Ro 
[Equations (62), (63), and (64) ], while the shape of the response curve 


depends on the ratio of (wM) 2 to one-half the sum of the squared 


resistances [Equation (65) ]. 


Example 


At this point it is of interest to look at some examples of 
coupling for a radio-frequency transformer. A variable transformer 
with a fixed primary circuit and adjustable secondary inductance and 
coupling coefficient will be used. Operation will be at one fre- 
quency, but the secondary load will be different for the various 
cases. The basic circuit is shown in Figure 22. Parameters for the 


transformer primary are given below: 


L, = 430 uhenry X, = j41.9 ohms 
1 Ly 

R) = 4,2 ohms QQ, = 10 

C, = 0.25 ufarad f = 15.5 kHz 


Table IV shows the transformer response for four distinct 
cases. In Case One the secondary load resistance is 20 ohms, in 
Case Two it is cut in half to 10 ohms. For Cases Three and Four the 
secondary load has been reduced .to one-third and one-fourth of its 
initial value respectively. 

In all cases the secondary Q is very low, indicating a wide 
bandwidth. Coupling coefficients are given as determined by Equations 
(59), (60), and (71). The latter equation is considered the most 
accurate, since it is based on measured rather than calculated values. 


Capacitor C, is used to tune the secondary alone to resonance at the 


2 


operating frequency. 
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Figure 23. Resultant Circuit 
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TABLE IV 


TRANSFORMER COUPLING ARRANGEMENTS 


Case 


Parameter 
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For Case One the coefficient of coupling is less than critical, 
Equation (63) shows slight undercoupling and the shape of the secondary 
response curve should be single-peaked according to Equation (65). 

Case Two is overcoupled by Equation (64) and will have a double-peaked 
curve from Equation (65). By the same equations Case Three is over- 
coupled but should have only one peak, a case of transitional coupling. 
Case Four is also a transitional form of coupling which should yield 
only one peak for the secondary curve. 

Using Equation (53) the transformer in each of the above cases 
may be replaced by a complex impedance. But since the secondary is 
tuned to the operating frequency only a pure resistance will be seen. 
The primary is also tuned, so the resultant load for the generator- 
resistance source is solely the reflected resistance, marked as R3 an 
the Table and in Figure 23. Only in Case One is the load less than 
the primary resistance; in the other cases there is a definite over- 


coupling effect. In summary, the only noticeable bandwidth limiting 


will occur in Case One. 
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Dr Signal Transmission at VLF 


As discussed above, VLF signals have predictable propagation 
characteristics, with a dominant ground-wave transmission mode. They 
may be received reliably at long distances by aircraft, surface sta- 
tions, and submerged vehicles. For these reasons present day use has 
been generally dominated by the military and scientific communities. 
High radiated power is the basic requirement for successful operation. 

The entire VLF band from 10 to 30 kHz is only 20,000 Hertz 
wide. It is important to use this resource in the best manner pos- 
sible, which means a tradeoff between the number of channels available 
and the bandwidth of each channel. Since the amount of information 
and the speed at which it is sent are determining factors in selecting 
the required bandwidth, it is important to consider all the factors 
involved. 

In the discussion of VLF aerial systems it was pointed out 
that the aerial itself is usually a compromise in terms of the physical 
size which can be built. Due to top-loading and the small radiation 
resistance normally encountered the bandwidth is often less than 100 
Hertz. (For the aerial discussed in the Impedance Matching section, 
the aerial bandwidth alone was about 38 Hz while the overall aerial 
system bandwidth was 76 Hz.) Although the overall system bandwidth 
may be increased (as discussed in the following section) there is 
usually a high price tag attached. 

Prior to the 1950's most signalling at VLF was by CW (continuous 
wave) emission, also known as radio telegraphy. Reception was accom- 
plished by human operators who transcribed the CW signals into letters 


and words. Because of this keying speeds were generally limited to the 
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range of 15 to 20 words-per-minute, and the small bandwidths required 
were well within the typical 100 Hz aerial system limitations. In 

fact, much effort was expended in trying to reduce the loss resistance 

in the system so that power could be more efficiently radiated. The 
corresponding reduction in. bandwidth. [see Equations (1l.a) and (2) ] 

would merely reduce the available system response down to that compatible 
with the keying requirements .* 

It is important to examine the more common form of emission 
used in radio communications, especially those which enjoy wide usage 
in the High Frequency (HF) spectrum from 2 to 30 MHz, and determine 
their applicability to VLF signal transmission. 

The normal voice or audio channel is considered to have a 
3,000 Hz bandwidth requirement. Standard double-sideband Amplitude 
Modulation (AM) has two identical sidebands, so each station using AM 
would require at least a 6,000 Hz.channel. This mode of emission may 
be eliminated from consideration since only three stations could 
operate in the entire VLF band, and. with. the aerial system discussed 
thus far the 100 Hertz bandwidth restriction would not permit 
operation. 

Single-sideband (SSB) with or without. carrier,. can’ be considered 
to have only a 3,000 Hz bandwidth requirement. Although the number of 
operable channels could now be increased to six, the aerial bandwidth 
restriction would still not allow proper operation. 

Narrow Band Frequency Modulation (NBFM) systems are usually 
considered to have roughly the same bandwidth as normal AM systems, 


when carrying an audio signal, while Wide Band Frequency Modulation (FM) 
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is by definition a system with a greater modulation index (and 
consequently greater bandwidth) than NBFM. Both of these systems 
may be eliminated from consideration. 

The above systems (AM, SSB, NBFM, and FM) have been discussed 
using an audio signal. In the HF bands, it is common practice to use 
frequency multiplexing and place a number of Teletype channels in the 
place of one audio channel. The data rate of such a system is much 
greater than the original voice signal, and makes good use of a given 
3,000 Hz bandwidth. It would be interesting to try this same approach 
at VLF, but once again the determining factor is the bandwidth of the 
aerial system. 

Perhaps the most satisfying results may be obtained by starting 
with the major limitation, the aerial bandwidth, and seeing what type 
of system can best be used and still meet the restrictions. If 100 Hz 
is taken as the bandwidth for each channel, as many as 200 channels 
may be put in the VLF spectrum. 

Since CW operation was the most common form of signalling prior 
to the 1950's, it is important to look. at its operating characteristics 
in more detail. It is basically a form of Amplitude Modulation, where 
the modulation is set at 100 per cent; that is, the carrier is either 
full on or completely off. This is binary (two-state) keying in its 
simplest form. 

Each letter and number is made up of from one to five parts or 
units. The basic unit is the "dot," a short duration "on" pulse. A 
"dash" is a longer pulse, equal in time duration to three consecutive 
dots. Spacing between dots and dashes is an "off" pulse (actually no 


pulse at all) equal in length to a dot. An "off" pulse equal in 
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length to a dash is used for spacing between letters. See Figure 24. 
The basic word has six characters, five letters and a space, plus two 
extra dots at the end of the ee The standard length word for 
determining signalling rates in the military is CODEZ, which is the 
equivalent of 60 dots long. 


If a word consisting of five consecutive letter "e's" was to 


" ! 


be sent it would require five dots for the "e's" alone and four dashes 
for the spaces between the letters, taking a total time equivalent to 
17 dots. The numeral zero (%) is represented by five dashes, with a 
space of a dot between each dash, for a total time of 19 dots. When 
sending a word consisting of five zeros, a time corresponding to 107 
dots would be required. 

In using CW consisting of only letters and numerals (punctua- 
tion marks are usually six units long), a single word may range from 
17 to 107 dots in length. The International Morse Code is arranged 
such that the more commonly used English letters are shorter in time 
than the less common letters, but even so, one word may be more than 
six times longer than another. 

The bandwidth requirements for CW using International Morse 


Code can be expressed approximately Fogg 


eee 
BW ow +s WPM cet ey 
where BW aw = bandwidth in Hertz 
WPM = Morse Code words-per-minute 


The standard word CODEZ has been used to determine Equation 
(72), so the expression is applicable to both English words and random 


"five letter code groups," since it represents an average word length 


of 60 dots. However, the "five letter code groups" (often used as a 
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[Letter"A" is dot-dash] 


Figure 24. International Morse Code 
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one letter "A" 


Figure 25. Teletype Character (Letter "A") 
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basic form of intelligence coding) will be random in terms of the 
letters used, and thus show more variation in length than a normal 
English text. The resulting wide range in sending time causes a loss 
of effectiveness in such a coding procedure. 

In the last ten years machines have been developed which will 
"copy" CW transmissions, replacing the normal human operator since they 
are ie a accepting speeds which before only the most tdlénted 
operator could copy, and with far fewer errors. Digital logic circuits 
are employed and once properly tuned the machine‘will produce an output 
printed on a normal Teletype page printer inserting carriage return and 
line-feed signals at the appropriate times. With a’ Boehme mechanical 
CW keyer at the transmitter CW rates at up to 60 words-per-minute may 


{ 


be used. 
a | Averena line method in widespread use today ‘ts Teletype, where 
Ene eee anata teen éaploved ‘are typéwriter-like keyboards’ and 
printers. ia mR reeettte each letter, tuméral and pdadr@ticn mark are 
of Seva dengeh, Wie that any Conbindtion of five letters ?-nimé@fdls or 


3.4 in 


Baa oan marks will tai exactly. the sdine timé’to send" The basic 


“sf A 


system uses a seven unit character. ‘The first’is” ‘a "stati" unit, the 


: fe > Daas 3 ho 2 Ass ‘ ‘= | é or "sc x | z 
next five depend on the element. to be sent, and the last‘ is’ a’“stop" 
IIo Ys 


unit. 


wh 
‘ i ' @ *, ae 
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Both Teletype. and CW. Morse. code. are basically binary’ signalling 


u ‘ 4 2 i ak 


— In the first, however, each character always has seven basic 


¢ 
t 
‘ 


units, net in telegraphy ie may ‘vary ‘from one to’six. ‘The’ five 


Ei earersieie units in rte must “alway besent, but Since they 


nS 


are aid of equal length, the arrangément of the five units-does’ not 


effect the time required to send them, as is the case"in CW: ‘The 


ilk 


"dots" and "dashes" of telegraph sending are "on" aignals, with “off" 
pulses used for separation. In Teletype the five information units are 
also distinguished by the difference of two states, but the states are 
not strictly. limited to "on" and "off." The reference level may be 
set at the bottom, as shown in Figure 25, giving on und off pulses 
(called "mark" and "space" respectively), or it may be set to the 
upper limit. Either of these systems, using positive or negative= 
giving pulses, is known as "polar" keying. I£ the middle of the pulse 
4s chosen for reference then the keying is called "bi~polar." 

With five "information" units, each of which may be in either 
of two states, there are 2? or 32 possible combinations. The 26 
letters of the English alphabet may be included, but there ts not room 
for all of the numerals of the decimal system. To satisfy this require= 
ment two additional characters are used: "Figures" and "Letters," The 
first shifts the keyboard arid printer to what on the normal typewriter 
would be the upper case position, Tlie characters used for letters in 
the lower case or "Letters" position may fiow represent tiumefala, punetua= 
tion marks, weather symbols, or special symbols such as a beli, if a 
“Letters” character is bhadn transmitted a shift is made to lower case 
and alphabetie ‘tencters &Y@ Once Hore sett. 

thus with seven-unit Teletype coding up to 62 different eharac= 
ters, each requiring the seme tiwe duration, may be sent, | 

In any binary eystem it 1s diiportant to be able to properly 
recognize which of the two states 16 belig received, With QW operation 
the two states usually chosen are on and off. Slice reception 16 fost 
ofter done by a2 huwan operator, the relative spacing between tha lengths 


of dots and dashes can be varied sevewhat and perfect copy st1ii reguit, 


Teletype operation, however, depends on the timing of the pulses, and 
the allowable variation from the proper timing is much less. Addition- 
ally, Teletype signals are decoded by machines, so the adaptability of 
a human operator is not available. 

If the binary levels of zero and.one are chosen, it is a simple 
matter to distinguish the two signals. When noise enters the system 
the basic reference level, if initially chosen at zero, will become 
instead the noise level. To counteract this a "threshold" may be 
chosen, say at 4s unit, so that any pulse received which is greater in 
magnitude that ¥ will be considered as a signal of one unit; signals 
less than 4 will be considered as zero level. Noise spikes of magnitude 
greater than 4 may cause an error,, but the average noise level will not. 
The proper choice of threshold level is not as simple as presented here, 
and quite a bit of work has been done on coding and levels, especially 
at viele 

Rather than. a. choice between two levels, improvement in the 
decision process could be realized if there were.more accessory require- 
ments. For example, a Teletype mark could be transmitted at one radio 
frequency and the space at a second frequency. Then a true mark would 
require a signal to appear at the proper mark frequency and no signal 
to appear at the space frequency. Detectors at the mark and space 
frequencies would employ the threshold logic described above. 

Table V shows the four possible mark-space combinations in a 
noise-free environment. If the transmitter is properly operating 
Cases One and Three would never be received. The threshold detectors 


might produce such an output, but sophisticated circuitry would remove 
them, or else properly indicate an error. Cases Two and Four would be 


proper signals. 
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TABLE V 


MARK-SPACE TRUTH TABLE 


Comments 


Transmitter off, severe signal 
fading. Logic circuit would 
recognize as incorrect signal. 


Space. 


Incorrect signals, logic would 
recognize and indicate as an error. 


Mark. 


Error correcting codes and similar plans to improve accuracy 


have been proposed and are in use. But the problem at hand is to make 
use of two separate signals for..Teletype signalling at VLF. 

(It is seen that CW could be used with this arrangement, sending 
separate dots and dashes. But the very attributes pe the Morse code, 
1.e., variation in character length, are a detriment here, especially 
in terms of selecting proper timing intervals.) 

Two methods of generating separate mark and. space signals for 
radio communications are commonly used. The first is Frequency Shift 
Keying (FSK), where the transmitter carrier is shifted between two 
discrete frequencies. Audio Frequency Shift Keying (AFSK) has the 
carrier transmitted continuously with two signals present in an audio 
sideband: one audio signal used for mark and the other for space. Where 
more than one Teletype channel is to be sent the AFSK method is generally 
employed, for here there may be several pairs of mark and space signals, 
all within a 3,000 Hz audio sideband. The carrier, in fact, need not be 


sent at all if a single sideband unit is used, and for elaborate systems 
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it is possible to transmit up to four independent sidebands from the 
same transmitter, each 3,000 Hertz wide, and each containing 16 or 
more separate Teletype channels. 

In the VLF spectrum phase-continuous shifting of the main 
carrier is the most advantageous method of employing rsx 7° 

Phase continuity is required to maintain a high level output 
and avoid serious transients in the antenna system, and in effect 
means that when shifting from one frequency to another both waveforms 
must be the same, i.e., both maximum or both passing through zero in 
the same direction. This may be achieved by using a single variable 
oscillator, switched at the signalling rate, or by switching the output 
between two separate oscillators. The difference between the mark and 
Space signals depends on the system employed; at HF frequencies it is 
commonly 850 Hertz. A radio system meeting the phase-continuity 
requirements is actually using a form of Frequency Modulation (FM). 

In simple terms, an FM system is initially at rest at some 
assigned frequency (carrier or center frequency)... If. a signal such 
as a sine wave is to be transmitted,. then the. radiated frequency will 
shift up from center frequency as the magnitude of the sine wave 
increases from zero, reaching its maximum deviation from the assigned 
carrier frequency when the sine wave is at its maximum value. As the 
magnitude of the sine wave decreases the carrier will be lowered in 
frequency until back at its assigned frequency. At this point the input 
sine wave will be at zero magnitude. The carrier will then be shifted 
down in frequency as the sine wave goes negative, finally returning to 
the’ center frequency at the end of one cycle. See Figure 26. The 


maximum deviation of the carrier from the center frequency is a function 


15 
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Figure 26. Simplified FM Frequency Plot 
Amplitude 


A. Equal time periods 
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B. Actual Teletype Signal 


Figure 27. Theoretical and Actual Teletype Waveforms 
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of both the magnitude of the sine wave and the modulation index of 
the FM system. 
The "modulation index'' is defined as: 


Af 
ra O73) 


Ss 


where Af = deviation of upper (lower) frequencies 
from the center frequency. 


BE 
Ss 


signalling frequency. 

The maximum deviation allowed by the transmitter is a specified 
parameter, so the modulation index normally varies with the information 
signal. By definition Narrow Band FM has an index much less than > 
while normal or wide band FM has an index equal to or greater than 5. 

For Teletype signals the carrier frequency is usually shifted 
upward for a mark and downward for a space, in both cases using the 
maximum allowed. deviation.. Since the shift. is not. instantaneous the 
signalling waveform may be "shaped" to reduce keying transients. 

In Teletype operations it is common to speak in terms of a 
data rate in words-per-minute. For military. usage, both 60 WPM and 
100 WPM are common, with recent emphasis on the latter. The equations 
below are based on the normal word consisting of six characters (of 
seven elements each), five letters and a space. The elements per 


second are of most interest, however, and for this the term "baud" is 


used. A baud is the basic data or element rate, that is, elements per 


eecends 
WEM = (elements ea character, 1 word. | 60 seconds, 
second 7 elements 6 characters minute 
WPM = (1.43) x Baud » pete7.4) 
Baud = (0.7 ) x WPM oo) 
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From Equation (75) 60 WPM Teletype has a baud rate of 42, so 
each identical pulse will have a time duration of 23.8 milliseconds, 
where time is the reciprocal of frequency. This is shown in Figure 
BIA). 

In practice the pulses are not of equal time duration. The 
stop pulse is made 1.42 times longer than the length of the start and 
information pulses, so the standard character requires 7.42 units of 
time. This means that each time unit must be approximately 22 milli- 
seconds long and the stop pulse will be approximately 31 milliseconds 
long. The average length is still 23.8 milliseconds, but the shortest 
pulse present is of the most concern and sets the basic bandwidth 
requirements. 

A signalling speed encountered where United States and allied 
military communications are concerned utilizes a 20 millisecond 7.42 
unit code. This corresponds to 50 baud operation, and by Equation (79) 
has a data rate of approximately 67 wem. +? 

Fourier analysis. shows that an ideal square wave consists of 
an infinite number of sinusoids at various magnitudes. In the FM FSK 
method of transmission the signalling frequency must be known to 
properly determine the modulation index and frequency deviation 
[Equation (73)]. As a first approximation a sinusoid may be fitted to 
the square wave which has the same magnitude, zero crossings, and 
period. The resulting waveform then has a period of 40 milliseconds 
or a frequency of 25 Hz, as shown in Figure 28. 

The signalling frequency of 25 Hz is a rough estimate at best, 


but with it an idea of bandwidth requirements may be ceeeeiaaea- 
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Figure 23, Pulse, Square-wave and Sine-wave 
As Signalling ‘aveforms 


79 


BW= 2 (f, +f) ees 6) 


where BW = total bandwidth required 


Fh 
1] 


frequency deviation for mark or space 

f = signalling frequency 
Proper solution of this equation requires knowledge of the allowed 
deviation. Studies of frequency deviation have been made through a 
study of the modulation index. Indices of one-half and one have been 
found particularly useful for square wave keying. An index of one-half 
uses all the signal power from the transmitter to carry information, 
but there is no discrete power at any one frequency and synchronization 
is difficult. When an index of one is used there are discrete spectral 
lines at the upper and lower shift frequencies, but they also take 
one-half of the signal power and carry no message information. ?° 

In actual message transmission random pulses rather than square 
wave keying will be useful. In this case the spectral density will 
vary from those described saaveie® Here the choice of a modulation 
index of one for phase continuous. binary. keying appears to be optimum. 

The application of FM to high power VLF transmitters has been 
successfully accomplished only within the last 20 ee Early 
studies in keying low frequency transmitters using reactance tubes were 
carried out at Westinghouse Corporation in early 1950. An experimental 
FSK system using a baud rate of 50 with a modulation index of one was 
operated between San Diego, California and Leesburg, Virginia with a 
radiated power of one kilowatt in the early 1960's.7+ 

One difficulty in using FSK at VLF is the small frequency 
shift (25 Hz for a modulation index of one for the 67 WPM discussion 


above) as opposed to the more common 850 Hz shift in the HF region. 
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Such small shifts can be successfully detected as reported in one 


study which used only 18 Hz for 100 WPM paeeoticaeeas 


E. Broadbanding 


Narrow bandwidth of VLF aerials is the major limitation to high 
speed information transmission. The impedance matching and transformer 
circuits have been shown to be practically lossless, and since the 
bandwidth (and Q) of the aerial is directly dependent upon the resistance 
in the circuit, an increase in resistance will increase the bandwidth. 
At the same time, however, there is a corresponding decrease in system 
efficiency. 

A direct attack on the problem would be to place a fixed, high 
wattage resistance across the parallel circuit, as shown in Figure 29. 
The new resistance will lower the circuit Q and increase losses. The 
variable inductance may be used to vary the influence of the resistor 
on the circuit, giving the effect of a variable resistance. 

In practice such a broadbanding resistor would be of comparatively 
high value and with the series inductor set at its full value have only 
a slight effect on the original antenna circuit. Since the antenna 
circuit in general would already have been designed, the best position 
would be between the aerial coupling unit (helix system) and transmis- 
sion line impedance matching network. Then the latter network could be 
slightly retuned to maintain a proper input impedance. See Appendix A. 

A second approach is more sophisticated and will not introduce 
as much loss (in theory no loss at all!) but calls for special devices. 
The idea is to automatically tune the aerial system to resonance in 
synchronism with the information signal. Under normal conditions the 


aerial has a very sharp response curve with 3 db points less than 50 Hertz 
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Figure 30. 


Switchable Reactance For Broadbanding 
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either side of the center frequency. If a bandwidth of 150 Hertz was 
desired, and FSK was used, then in theory the antenna could first be 
tuned for operation 75 Hz below the center frequency; then with a new 
reactance in the circuit be retuned to resonance 75 Hz above ‘the center 
frequency. If the new reactance was switched at the same rate as the 
FSK, then at both mark and space frequencies the aerial would be 
resonant. 

Figure 30 shows this schematically, using a reactance jX in 
parallel with the tuned circuit. At the lower frequency fi the basic 
circuit components are used; at the high end (f,) the new reactance is 
added and it is tuned for circuit resonance at this frequency. 

In 1954 Jacob and Bratich put forth the idea of using a saturable 
reactor for the switched tuning Semone oo An experimental model of 
such a device was found to dissipate about 20 KW when full frequency 
shift was employed using a 500 KW transmitter. Special modulation 
equipment was required to obtain the desired results, and the saturable 
reactor itself measured over 1,400 cubic feet and weighed about 23 
tons. Operation with Teletype FSK was successful using a 100 Hertz 
shift at 19 kHz. 

Wolff reported in 1957 on a variation of the method above, 
using a variable capacitance as the switched reactance. Vacuum tubes 
were keyed with the Teletype signal and effectively varied the 
reactance. This method was applied to a 10 KW transmitter operating 
at 28.5 kHz using a shift of 100 Hz. However, modulator tubes for 
high power (above 10 KW) VLF operation were not sven ie’ 

More recently the problem of transients in the tuning arrange- 


ments above have been considered. Galejs points out that in an LC 
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resonant circuit the energy is alternately stored in the inductive 
and then the capacitive elements. If one of the elements (or part 
of that element) is switched out of the circuit at the instant it 
contains all the energy available, there will be no transients. The 
element may then be switched back into the circuit when the other 
reactance has again reached a zero energy level. Although the author 
discusses FSK applications at VLF and variation in parameter values, 
switching timing and the switches themselves, there appear to be many 
practical problems which would require a great deal of Seidvinaed za 
The ideas expressed by Wolff and Jacob generated a great deal 
of thought in the area of reactance variation for increasing bandwidths 
in linear circuits. The use of principles normally associated with 
servomechanisms and magnetic devices were applied to the problem by 
Johannessen in Mae He suggested the use of two non-linear 
inductors connected in series with the VLF aerial tuning coil. The 
average inductance of these devices is caused to vary in accordance 
with the keying signal, through the use of a feedback signal derived 
from the phase difference between the antenna current and transmitter 
voltage. No computations regarding the size of the non-linear induc- 
tors are given, but additional equipment for phase detection and 
inductor control are required. The author states that a small-signal 
transient analysis shows the system to be limited to a new total band- 
width equal to one-fourth the bandwidth of the power amplifier which 
drives the inductors. This, however, should not pose any real problems 
at the VLF frequencies where the system bandwidths are relatively small. 
The most promising idea for increasing aerial bandwidth was 


presented by Zverev and Blinchikoff in teeuees Here the aerial radiation 
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resistance is considered as a resistive load and the transmitter as a 
source with series resistance. The desired bandwidth is determined as 
part of an overall frequency response and by network synthesis a cir- 
cuit having the desired Butterworth or Chebyshev response is designed. 
The aerial reactance is then included in the analysis and combined with 
additional components to complete the bandpass filter network. The 
article gives numeric tables for circuit design including response 
time, ripple and bandwidth. When a two-pole Chebyshev filter is de- 
signed (which is recommended as the most economical) the resulting two 
resonant peaks may be adjusted to coincide with the mark and space 
frequencies for FSK Teletype transmission. 

The example filter included by the authors was for pulse or 
Square wave modulation and was built as a simple Tee network. Although 
impedance transformation was not a primary concern in this case, it 
could perhaps be included in a more thorough analysis, or else a 
separate wideband impedance matching network could be added. 

When a VLF aerial is connected to a. transmitter an additional 
resistance will be placed in the circuit. since the power amplifier may 
be represented as a zero-impedance source in sesetlee with its internal 
resistance. Thus there will be an increase in bandwidth. Since power 
amplifiers for VLF are usually operated Class C, the exact value for 
the generator resistance is difficult to determine. It has been sug- 
gested that as a rough guideline the bandwidth of the antenna system 
can be increased by an amount found by multiplying the present band- 
width by the reciprocal of the Class C amplifier efficiency. If this 
procedure proves correct the tubes could be operated at lower effi- 


ciencies and the increased bandwidth achieved. This method is 
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expensive in terms of power dissipation and overall system efficiency, 
but an increase of approximately 12% can be achieved if the amplifiers 
operate at their nominal 80% efficiency. 

In summary, some form of broadbanding for the normal VLF 
antenna system is required if high speed (wide bandwidth) information 
transmission is desired. This may be done by brute force -- the 
introduction of additional resistance into the circuit. Less power 
will be lost if a variable reactance is switched into the aerial cir- 
cuit synchronized to the signalling information to tune the aerial to 
the separate mark and space frequencies. Transients are a problem 
here. The most potentially beneficial broadbanding method is to 
synthesize a bandpass filter having the desired frequency response. 
Regardless of which method is used, an increase in bandwidth of the 
overall system (but not the aerial itself) will be provided by the 


source resistance of the power amplifier tubes. 
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bis SYSTEM USED AT NORTHWEST CAPE, AUSTRALIA 

In the Basic Principles section VLF aerials, matching and 
coupling networks, signalling methods and bandwidth widening tech- 
niques were discussed. For any particular VLF installation all of 
these areas will be important, but certain areas will vary among 
installations. In the writing of a computer program for VLF systems 
numerical values must be used, at least in the checkout phase of the 
program. 

Rather extensive data on existing U. S. Navy VLF stations, 
as of 1963, has been omupuiled However, the installation at 
Northwest Cape, Australia was under actual construction during the 
time this thesis was being written, and could be expected to incor- 
porate the latest techniques in the VLF art. Additionally, William 
E. Norris, Professor of Electrical. Engineering, Naval Postgraduate 
School, was assigned to the installation and could provide up-to-date 
information. Field engineers were making minor corrections and 
numerous measurements daily, so a great quantity of detailed informa- 
tion was available. For these reasons, the Northwest. Cape installation 
was chosen as the model for the computer program. 

The transmitter-antenna system was designed to operate in the 
frequency range from 15.5 kHz to 30 kHz, with a radiated power of 
1,000 KW. Major signalling is by 50 baud Teletype, using FSK with a 
shift of 25 Hz and a modulation index of one. CW operation. would be 
available on an emergency basis. The assigned operating frequency is 
15.5 kHz 

Table VI shows the major system parameters of the installation, 


and Figures 32 through 37 show the general schematics of the various 


units. 
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TABLE VI 


DESIGN PARAMETERS FOR NORTHWEST CAPE, AUSTRALIA 


General 


Aerial 


Helix 


Bandwidth 


Resistor 


Operating Frequency 


Radiated Power 


ILS) .,5 Uesle 


1,000 KW 


Vertical monopole, with top-hat loading 


Effective height 
Required efficiency 
Static Capacitance 
Downlead inductance 
Radiation resistance 
Loss resistance 
Self-resonance frequency 
Aerial bandwidth 
Bandwidth (with loss) 


Antenna current (approx.) 


Tuning Variometer 


Coupling Variometer 


Input impedance (at resonance) 


Resistor 
Maximum dissipation 


Series Variometer tuning 
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628 feet 
50% 

0.164 wfd 
139.0" un 
0.154 ohms 
0.154 ohms 
Jano KHz 
36.125 HZ 
76.250 Hz 


2550 amps 


484 uh 
20 uh 


12.5 ohms 


33.25 ohms 
500 KW 


200 - 800 uh 


TABLE VI (continued) 


Tee Input impedance 20.0 ohms 
Output impedance 12.5 ohms 
Leg One Inductor 259— 275 wh 
Leg Two Capacitor "06 -¥%.6 wéd 
Leg Three Inductor 50 - 300 uh 
Transformers | Output impedance 20 ohms 
Primary winding 430 uh 
| Primary capacitor OZ 5m ted 


Secondary inductance and capacitor settings vary 
with number of units being used. 


i One unit Inductance 73° uh 
Capacitance 1.475 ufd 
k 0.522 
Two units Inductance 49.0 uh 
Capacitance 1.076 wid 
k 0.474 
| Three units Inductance 35.3.0 
Capacitance 0.996 uid 
k 0.444 
Four units Inductance 27-0 in 
Capacitance Oi 3i/ Omilixel 
k 0.508 
Power Four each AN/FRT-67 power amplifiers 
Amplifiers rated at (each) 500 KW 
Desired plate-to-plate load 420 ohms 
Signalling Teletype 50 baud 
FSK, modulation index 1.0 
Required minimum bandwidth 100 Hz 


CW for emergency use 
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Figure 32. Aerial System 
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Figure 34. Bandwidth Resistor Network 
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EVs COMPUTER PROGRAM 
A. Outline 

A program written for a high speed digital computer can be 
made to carry out simple tasks in an extremely short time. By com- 
bining a number of simple steps a more complicated operation may be 
performed. When using a specialized program for aid in designing a 
VLF antenna system some basic assumptions are made and the various 
parameter values may be permuted to observe the effect on the system. 

The computer program written for this thesis (VLFANT) used 
the installation at Northwest Cape, Australia for the basic model. 
Thus the output or "results" apply only to this particular facility. 
However, VLFANT is sufficiently general that component values may be 
rapidly changed and the network may be modified with relative ease. 
In fact, a number of existing antenna systems can be run with this 
program by merely modifying aerial parameter values. 

VLFANT is written in FORTRAN IV language for use on an IBM 
System/360 high speed digital computer. Full use is made of complex 
variables. The circuit under study. is broken into five major sections: 
the aerial, the helix matching circuit, the bandwidth resistor circuit, 
the Tee matching circuit, and the power amplifiers-transformers 
coupling circuits. Each of the major circuits is put into a separate 
subprogram so that modifications may be made without effecting the 
remainder of the program. 

In addition to the five Major Subprograms above, the Main 
program acts as the control section. It "calls" the Major and Minor 
Subprograms (discussed on the next page) and by making simple changes 


in Main the program output may be varied to give the designer precisely 
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what he wants each time the program is run. Also, the entire program 
except for Main may be put into "object Deck Format'' which means the 
computer cards are changed from their original FORTRAN IV language 
into machine language. Then when the complete program is run only 

the Main program needs to be "compiled" (changed into machine language) 
yielding an important savings in computer time. 

Six Minor Subprograms are included in VLFANT. These sections 
operate on the major circuit sections by putting them into the network 
analysis "ABCD" matrices, determining network response, variation in 
circuit Q, currents, voltages, power and phase relationships throughout 
the networks as a function of frequency. This data may then be 
displayed in a number of ways, including graphically from the standard 
paper printer unit. Detailed descriptions of each of these sections 
are given below. 

Seven Accessory Subprograms are included. They are used by 
the subprograms above to perform repetitive mathematical manipula- 
tions, supply tuning limits for variable components, perform operations 
on matrices, and plot output data generated by other subprograms. They 
are described in detail later in this section. 

For flexibility variable names in VLFANT were chosen to match 
normal circuit theory as closely as possible within the limits of the 
FORTRAN IV language. A great emphasis was placed on indexing these 
variables so that new circuits could be substituted without creating 
numerous program corrections. Figure 38 shows with block diagrams the 
basic designation of the antenna system variables. Note that each of 
the five major circuits is numbered so that all measurements and para- 


meter values dealing with a specific section will have an index number 
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Aerial 
Network 


#5 


Power Tee Bandwidth Helix 
Amplifiers Network Network 
-Transforme #2 #4 


#1 


Figure 38. Major Circuit Sections 


COlmnGS. 1.) CAP(5,1) 
VOLTS(5,1) —— RAD 
AMPS (5,1) POWER(5, 2) 
POWER(5,1) PHASE(5, 2) 


PHASE(5, 2) 


Figure 39. Aerial Network 


oF, 


corresponding to that block. Each block is further divided into an 
input and oueut designation. 

As an example of the indexing, Figure 39 shows the Aerial cir- 
cuit with components and circuit parameters indicated as they are used 
in VLFANT. Since this block number is five, all variables have a "5" 
as the first index. Currents, voltages, and so forth dealing with the 
input. have a "1" as the second index; at the output (in this case the 
radiation resistance) the index is "2." Reactive elements are 
numbered sequentially, and reactance values follow the same notation. 
[XL(5,1) is the complex impedance of COIL(5,1) at the designated 
frequency of operation, i.e., jw. | Resistances are the only excep- 
tion to the indexing terminology; since there are only three in the 
entire antenna system they are given descriptive names. 

A certain amount of duplication occurs with this indexing 
system. At the connection between. the Tee network and Bandwidth 
resistor blocks, for example, VOLTS(2,2) and VOLTS(3,1) are equal, 
and so are AMPS(2,2) and AMPS(3,1). Flexibility is maintained, 
since in one study it may be of interest to allow different power 
levels to be radiated by the aerial and determine the voltage and 
current levels required at the bandwidth resistor, in which case 
the Tee network is not of immediate concern. Or the effect of the 
Tee network may be under study assuming a fixed voltage at its input 
and a constant resistance load, in which case none of the other cir- 
cuit blocks are required. 

Further. network manipulation is possible since the bandwidth 
resistor may be removed and a transmission line inserted, or the Tee 


network might be replaced by a Pi network. Providing the same notation 


is followed, these changes may be readily accomplished. 


98 


A great deal of use was made of "COMMON" statements throughout 
the program. These statements allow a block of storage in the memory 
cells to be used in common by many subprograms, substantially reducing 
the overall storage space required. Many terms are common to the 
majority of subprograms, such as VOLTS(a,b) and ZTOTAL(a,b), and 
via Common they are always available when needed. Additionally, 
variables in Common do not need to be handled through a calling state- 
ment when using a particular subprogram, so simplicity in use is 
achieved. 

General operation of VLFANT is divided into two major areas: 
the first based on an impedance analysis and the second on a voltage- 
current analysis. 

Beginning with the given aerial component values and radiated 
power, the entire system is analyzed to find the impedance values at 
the input to each of the sections at the specified operating frequency. 
Impedance values for the helix and Tee networks are already specified, 
so the inductance and. capacitance. settings required are determined. 

The plate-to-plate resistances for one, two, three, and four power 
amplifiers are the end result of this part of the program. 

When determining values for variable elements a notice is 
given if the required setting is outside the tuning range. Computation 
continues with the computed value, however, so that later modifications 
of the tuning range may be made if desired. 

With values for all components known, the program then uses 
the radiated power and radiation resistance to determine the "load" 
voltage and current. Through the use of "ABCD" matrices the corres- 


ponding voltages, currents, power and phase relations are determined 
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for each network working from the aerial to the power amplifiers. 
Impedances using Ohm's Law are also computed, as a check on the values 
previously determined. 

At this point the frequency may be varied symmetrically on 
either side of the operating frequency by using subprogram DELTA. 
Impedance variations are determined for each circuit block, as are 
voltage, current and power levels, based on the plate voltage for the 
power amplifiers determined above being considered constant. The 
above data is presented in tabular form; it may also be graphed by 


calling the appropriate subprogram. 
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B. Major Subprograms 


This is the major control section for VLFANT. All components 
are generally defined here: values for all inductors and capacitors, 
resistances, radiated power, desired impedance levels (such as at input 
to helix, at both input and output of the impedance matching network, 
power amplifier plate-to-plate, and bandwidth resistor), tuning 
ranges for all variable. components, breakdown voltage levels, and so 
forth. Additionally all subprograms are called from MAIN, so just how 
the program is used can be determined here. 

Other parameters are also defined, such as the operating fre- 
quency, the number of steps of each side of center to be taken when the 
frequency response is run, the size of each step, whether or not the 
bandwidth resistor is to be used in the analysis, and so forth. 

Once the overall program.is set.up for a particular antenna 
system many parameters will have been roughly determined. For example, 
the tuning ranges for the reactances in the Tee network can be determined 
by hand, while the precise setting will be done by the computer. Then a 
simplication is possible by putting these "constants" into the BLOCK 
DATA program described below. This leaves MAIN. free to carry out its 
main function: that of controlling the overall program. 


An example of one MAIN program is included in Appendix B. 


PRINT: 

This subprogram supplies the first part of the computer program 
output. It prints the designated center frequency and desired radiated 
power. For the aerial the predetermined radiation and. loss resistances 


and. top-hat capacitance and down-lead inductance are presented. In the 
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bandwidth resistance circuit the value and maximum power dissipation 

of the resistor are printed as well as the tuning range of the variable 
inductor. Minimum and maximum values for the coils and capacitors in 
the Tee network are then printed. For the toroidal transformers the 
primary inductance is given, and then the coupling coefficient and 
secondary inductance of each transformer when used alone or as a part 


of a two, three, or four unit combination. 


AERIAL: 

This subprogram is concerned only with the aerial system, with 
the radiation resistance considered as the output parameter or load. 
The down-lead inductance and top-hat capacitance are the reactive 
elements, the radiation and loss resistances completing the circuit. 
All values are given from the initial conditions specified in the main 
program. The aerial self-resonance frequency, the '"'Q"' and bandwidth 
(using only the radiation resistance) and the efficiency are computed 
and printed out. With a desired radiated: power (from the main program) 
the required current is found as are the voltage across the radiation 
resistance and the "top-hat" and "base" voltages, those which appear 
across the reactive elements. The voltage across and current through 
the radiation resistance, VOLTS(5,2) and AMPS(5,2) are used as 
references for later voltage-current manipulations. The "output"™ of 
this program is the impedance seen at the terminals, ZTOTAL(5,1). 


See Figure 39. for the circuit and terms used in this subprogram. 


HELIX: 
This subprogram uses the impedance of the aerial system as its 


load and ZIN as the desired impedance at its output terminals. The 
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program may use either of two methods to determine the values for the 
tuning and series variometers. 

Method One considers the entire helix coil as the inductance 
used to resonate the aerial circuit [cancelling out the reactive part 
of ZTOTAL(5,1)] and then finds. the tap point to give a resistive load 
of ZIN. 

Method Two uses an analytical expression where the coil is 
tapped so that the input impedance will be purely resistive with a 
magnitude equal to ZIN. Either method.may be selected from the main 
program; the values for the series and tuning variometers are given in 
printed form, and the issue ZTOTAL(4,1). 


See Figure 40 for the circuit and. terms. used in this subprogram. 


TEE: 

Values for. the. legs of. the. Tee.network are chosen using the 
principle of a quarter wave transmission line, that is, the square 
root of the product of ZIN and ZOUT. The legs are numbered in ac- 
cordance with Figure 35. Settings for the legs are easily determined. 

Following determination of component values the tuning ranges 
for each component is checked to see if any out-of-range values have 
been selected, and, if so, the printout indicates. whether the tuning 
is high or low and the value of the limit which has been exceeded. 

Next the completed network is connected with the helix-aerial 
combination as a load, and the input impedance ZTOTAL(2,1) computed. 


See Figure 41 for the circuit and terms used in this subprogram. 
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COIL(4, 2) 


ZOUT 

VOLTS(4,1) VOLTS (4,2) 

AMPS (4,1) AMPS (4,2) 
POWER(4,1) POWER(4, 2) 
PHASE(4, 1) PHASE(4,2) 

ZIOTALA 4.1.) COIL(4,1) ZTOTAL(A. 2.) 

Oo 

Figure 40. Helix Network 
COLL (23) COUR 2,2) 

Oo 
ZIN ZONE ZQUT 

VOLTS @251) VOLTS (2,2) 

AMPS (2,1) AMPS (2,2) 
POWER(2,1) POWER(2, 2) 
PHASE(2,1) ZTHREE CAP (2,3) PHASE (2,2) 
ZEOTAILC2. 1) ZTOTAL(2, 2) 


Figure 41. Tee Network 


104 


TANK: 

This subprogram computes the plate-to-plate impedance for the 
power amplifier tubes using one, two, three, and four separate power 
amplifiers. 

For a single power amplifier. the. impedance seen at the input 
to the Tee network is used as the load. Combined with the transformer, 
a single complex series impedance is computed. This is transformed to 
a parallel impedance, and the primary capacitive reactance added. The 
resulting impedance is that seen plate-to-plate. 

When two or more amplifier units are used, it is necessary to 
determine the effect of the modified load as seen by each unit. Using 
the procedure discussed in the theory section, the. primary impedance of 
each transformer is found for the various. combinations. As above when 
using only one unit, a value for the plate-to-plate impedance is 
determined. 

The printout shows the inductance value of the transformer 
secondary and mutual inductances, as well as equivalent primary 
impedance and resulting plate-to-plate load for use of one, two, 
three, and four power amplifiers. 


See Figure 42 for the circuit and terms used in this subprogram. 


RESIST: 

This subprogram is used to determine the effect on the "Q" 
of the antenna system of the resonant frequency as the coupling coil 
is set at its high and low values. The resulting minimum and maximum 
"Q" are then printed out. Provisions are made to use the subprogram 


with an arbitrarily selected value for the coil; the resultant "Q" is 
computed and sent.to the main program. 


See Figure 43 for the circuit and terms used in this subprogram. 
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Figure 42. Power Amplifier/Transformer Network 


COIL(3,1) 


RBW 


Figure 43. Bandwidth Resistor Network 
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C. Minor Subprograms 


MATRIX: 

Four terminal networks may be described by using complex 
two-by-two "ABCD" matrices, where the input voltage and current are 
found by multiplying the output terminal voltage and current by the 
ABCD matrix. If several networks are cascaded, then the input voltage 
and current parameters of the overall network may be found by multi- 
plying the load voltage and current by the individual network ABCD 
matrices in sequence. If a particular network is removed from the 
overall circuit, then by merely leaving out its matrix the circuit 
may be re-analyzed. 

In the Matrix subprogram each section of the antenna network 
is described in terms of its ABCD matrix. Individual matrices are 
written for one, two, three, and four power amplifier transformers. 
Values for all components are either given initially or computed in 
the subprograms. described above. 

Two matrices. are developed for each section of the network if 
required, so that either set of terminals may be used for the input or 
output. In the aerial matrix, the radiation resistance is considered 
as the "load" or output impedance. 

Frequency is the major variable in this subprogram and each 
matrix can be labelled and printed out if desired; the results are 


always returned to the main program. 


DELTA: 
This subprogram is used to determine the impedance of the 


various parts of the antenna system as a function of frequency. Com- 


ponent values are those computed at the center frequency. Deviations 
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from the center frequency are made equally both above and below; the 
number of steps (up to 240) and the increments of each step are 
determined by the user in the main program. Impedances at the aerial, 
helix, bandwidth resistor and Tee network are found at each frequency 
and put into arrays. Either of the two configurations for the helix 
matching systems may be used at the will of the user. Although values 
for the bandwidth resistor connected to the helix are computed, the 
user may choose to compute the Tee network response by considering the 
resistor as removed from the circuit. 

In the TANK subprogram values were found for the secondary 
capacitor which would give a pure resistance load to the power trans- 
former at resonance, In this subprogram these predetermined settings 
are used when two, three, or four units are put on the line. Then the 
plate-to-plate impedance is determined as a function of frequency when 
any combination of power amplifiers. is used. These values are also put 
into array form in the subprogram. 

An additional impedance array is computed and sent to the main 
program, but is not available for printed output in the subprogram. 
This is the sum of the aerial impedance and the reactance of the series 
variometer of the helix circuit. Defined as ZTOWR in the program, it 
will be used later to determine "antenna" current. 

Output from this subprogram may be used in several ways. The 
impedance arrays, frequency array, and normallized frequency array 
(where the center frequency is set equal to zero) are all sent to the 
main program for further use, The printout of results is in tabular 
form: frequency and the network impedances in one table, and frequency 


and the transformer results in a second table. These printed results 
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may be omitted entirely, both used, or either one used at the 


discretion of the user. 


QVARY: 

This subprogram. determines the variation of the aerial as a 
function of frequency over the VLF band from 10 to 30 kHz, in 500 Hz 
increments. First the Q and bandwidth (in Hz) for the aerial alone 
are computed. Then the loss resistance is included and the Q and 
bandwidth again computed under this "loaded" condition. The frequency 
and two pairs of computations are put into array form and returned to 
the main program. 

It should be noted that the radiation resistance and reactive 
elements of the aerial used in. this subprogram are those used in the 
main program. They are generally computed or measured. at the frequency 
where operation is contemplated, and will vary with frequency, but are 
assumed constant here. for purposes of examination over the VLF band. 

The results may be printed in tabular. form or surpressed if 


so desired by the user. 


ANTAMP: 

This subprogram will determine the change in various currents 
in the antenna system as a function of frequency. For use here, the 
voltage on the plates of the power amplifier tubes when full rated power 
is being radiated is used as a reference. (The particular voltage 
selected will depend upon whether one, two, three, or four units are on 
the line.) It is assumed here that the voltage will remain constant 
regardless of the frequency, a condition which is not in fact true, but 


may be logically made. (See the theory section under Transformer Coupling.) 
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The frequency is allowed to deviate equally both above and 
below the center frequency; the number of steps and the increments 
are those chosen when subprogram DELTA was called. 

At each frequency, the current drawn from the power amplifier 
tubes is determined by dividing the plate voltage by the transformer 
impedance determined in the DELTA subprogram array for the corres- 
ponding frequency. The MATRIX subprogram now is used, at the chosen 
frequency, to compute the voltage at the output terminals of the 
transformers; this is divided by the impedance presented by the Tee 
network. At the helix network the voltage at its. input terminals 
is divided by the sum of. the aerial impedance and the reactance of 
the series variometer (ZTOWR - discussed above). This is defined as 
the "antenna" current. At the radiation resistance in the aerial 
circuit both the current flowing through it. and the power across it 
are computed. 

Thus this subprogram, by assuming. a constant voltage on the 
plates of the power amplifiers recompute the voltage and currents at 
each frequency, as well as the output power. 

The results, in array form, are sent to the MAIN program, and 
may be printed out in tabular form. Two tables may be printed, omitted, 
or either selected for output. The first table includes "antenna," 
helix, Tee, bandwidth resistor current and frequency; the second has 
transformer, aerial, and radiation resistance currents, as well as 
radiated power and frequency. 

GRAPH: 
This subprogram is used to graphically display data which has 


been previously computed, and possibly presented in tabular form. For 
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example, it may be of interest to know the variation of the helix and 
Tee impedance as a function of frequency. Since these values are com- 
plex arrays, new arrays consisting of the real and imaginary parts are 
generated. Then this information may be displayed graphically by use of 
PTPLOT, using a separate plot for each parameter... The independent 
variable is the normallized frequency array since deviation from the 
center frequency is of prime interest. 

If desired, two or more arrays may be combined into a single 
plot, such as the real and reactive portions of a complex impedance, 
and presented using PLOTIT. This procedure presents both parts of the 
impedance on the same scale, which aids in determining relative changes 
as frequency is varied. 

ABCD: 

This subprogram makes. use. of the. ABCD matrices. for the different 
networks in the antenna system. The desired radiated power used in the 
AERIAL subprogram to determine the voltage across and the current through 
the radiation resistance, [defined by the reference VOLTS(5,2) and AMPS 
(5,2)], is the basic reference. 

Subprogram MATRIX is called using the center frequency, and each 
of the voltages and currents for the various. networks. is found by appro- 
priate multiplication of matrices. Then subprogram COMP is used to de- 
termine the impedance, phase angle, and power at. each network, down to 
the input to the Tee network. This data is then printed in sequence 
starting at the radiation resistance. 

Next thé four different transformer combinations are connected 
to the Tee network. The parameters of voltage, current, impedance, 


phase angle, and power are determined both at the input to the 
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transformers and with the primary tuning capacitor connected, for all 


four combinations. 


D. Accessory Programs 


These programs are incidental to the main program in that they 
are used for setting initial conditions, doing special arithmetic manip- 
ulations, or displaying certain data. No direct. printed output is 
available from these programs unless specifically noted; the data is 


used in intermediate steps of various subprograms. 


BLOCK DATA 

This is a special program used to initialize many secondary 
variables. In particular, the tuning limits. for. variable inductors 
and capacitors are inserted here’ since. they will be "constants" for 


a given system. 


FACT1(A) 
This subprogram computes the sum of A squared plus one, 
which often arises in converting from series..to parallel circuits 


and vice versa. 


FACT2(A) 
This subprogram computes the quotient of A squared plus one 
quantity divided by A squared. This is also used in series-parallel 


transformations. 


COMP 
This subprogram is called with a complex. voltage and current. 
If either is zero a statement is printed so indicating; otherwise, the 


complex power and impedance are computed, as is the impedance phase 


angle, using current as reference. 


LZ 


MULT 

This subprogram is called with complex matrices (arrays), up 
to order 9 by 9. If the orders are acceptable for multiplication 
(a notice is printed if they are not), the resulting product matrix 


is sent back to the main program. 


PIPLOT 

This subprogram is a platting routine using the computer page 
printer and paper for output. With two given, equally dimensioned 
arrays of real numbers, one is plotted against the second, using 90 
horizontal spaces on the output paper. The values of the two arrays 
corresponding to each plotted point are printed. The two extreme 
values of the dependent variable are printed, and. titling for the 
column listing of the variables (such as "Frequency" for the 
independent and "Response" for the dependent) is available. If 
either array passes through zero value an. appropriate axis will be 
printed. The length of the printout (number of independent points 


plotted) is limited solely by the storage size of the computer. 


PLOTIT 

This subprogram will plot up to nine dependent arrays of real 
numbers against one independent array. Basically an extended form of 
PTPLOT above, this routine uses a horizontal field of 100 points on 
the page printer, since only the independent variable is specifically 
listed. A common scale for the dependent variables is used (allowing 
comparison of various responses) with values marked at the beginning 
and at 25% locations across the bottom of the plot. Axes, like those 


in PTPLOT, are available. The independent array may be sorted into 


ascending order if desired. Titling is available, such that eacn curve 
may be labelled along with its "mark" on the plot. The total number 

of independent variable points which may be plotted depends on the 
number of dependent arrays used. All arrays (including the independent) 
are formed into a single array within the program, which is limited to 
1,000 points. Thus three dependent arrays of length 250 may be used, 


or nine dependent arrays, each of size 100. 


E. Variables Used in VLFANT 

The choice of variables for VLFANT was made with quick recogni- 
tion as the main objective. Since names are limited to six alpha- 
numeric characters, in some cases a compromise was necessary. (This is 
the reason "aerial" is used so much in place of the more common 
"antenna"). 

All variables are single precision (16 bits) rather than 
double precision (32 bits). Although there is some loss in accuracy 
since numbers of quite different magnitudes are multiplied together 
(for example inductive values with orders of magnitude of iio” are 
multiplied by the radian frequency with order of magnitude of 10°) 
the loss is in the least significant digits and is not felt to warrant 
the use of double precision. 

Impedances are all complex, with the exception of pure 


resistances when they are "given" 


variables. ‘The complex "s," 
from Laplace analysis, is used extensively. 
Table VII lists the most important variables, Of particuler 


interest are those which must be supplied by the program user, since 


it is here that changes in the system must be made. <A number of 


i14 


variables not listed will be found in VLFANT, but these are generally 
used for intermediate steps and are of major concern only when making 


a detailed study of the programs. 


a 


TABLE VII 


VARIABLES USED IN VLFANT 


Defined By Units Comments 
U - User 


cycles } Center (operating) 


frequency 
WATTS R U watts Desired radiated power 
j at center frequency 
RAD R . U ohms Aerial radiation 
resistance 
| RAS R U } ohms Total aerial resistance 
(radiation plus loss) 
| CAP(5,1) R U | farads | Aerial static 
) | | capacitance 
COILGa) R U | henrys { Aerial inductance 
ZOUT R | U | ohms | Helix input impedance 
at operating frequency 
RBW R | U : ohms Bandwidth resistor 
; : | resistance 
} q 
PRBW R U | watts Bandwidth resistor 
; | dissipation 
ZIN R ; U ‘ ohms Tee input impedance 
at operating frequency 
RBWCH R i U { henrys| Bandwidth system vario- 
meter -— upper limit 
RBWCL R ; U i henrys} Bandwidth system vario- 
; meter - lower limit 
CP21H R U farads Cap(1,2) 
i Transformer Jupper limit 
CP21L R U farads| 9" °once" | Wigan Cas 
ower limit 
CP23H R U | farads} Tee Leg two ,Cap(2,3) 
upper limit 
GP23L R | U farads Cap(2,3) 


lower limit 
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eee pees 


CL21H R U henrys |Tee Leg Coil(2,1) 
one upper limit 


CL21L R U henrys Conveze 
lower limit 
CL22H R U henrys Leg Coil(2,2) 
two upper limit 
Gie2.2]. R U | henrys Gonl(2.2) 
lower limit 
GOTL C2ygls) R U henrys |Tee Leg one 
Value for coil in Tee 
section 
COTL C11) R U henrys |Inductance for primary 


of each transformer 


CLP (a) R U 


block, "b" input®or 
output) 


- Transformer coupling 
coefficient (a goes 
from 1-4) 
} TAP (a) R U | henrys }Transformer secondary 
lwinding (a goes from 
| 1-4) 
WRPTP ae U ohms Power amplifier plate- 
to-plate resistance 
ff CAP(1,1) : R U farads |Transformer primary 
i circuit capacitance 
COIL (351) Pa U henrys {Bandwidth resistor 
| coupling variometer 
i K R U - Number of steps in fre- 
| ;quency deviation program 
DEL R U - ,Size of each step in 
| [frequency deviation 
5 
OMEGA . P radians/|Frequency in radians 
‘ | | sec. | 
Is C Pp radians/|Complex operator, fre- 
sec. quency in radians 
VOLTS (a,b) c 3 volts Voltages ("a" gives 
block, "b" input*or 
‘ output) 
| a 
AMPS (a,b) C | E amperes |Currents ("a"' gives 
: : 


nN 


biliock; sb input or 


ee eee 


j 

ZIOTAL Ga, b) P ohms Impedance ("a" gives 

output) 

YTOTAL(a,b) C P mhos Admittances ("a gives 
block, (by inputior 
output) 

POWER(a,b) C P watts Power ("a"' gives block, 
"b'' input or output) 

PHASE(a,b) C | P degrees! Phase angle ("'a" gives ‘ 
j (rad.) ieee’ 'b” input ior 
i ' output) i 
y 1 f 

COIL(a,b) i R } P henrys {Inductors ("a"' gives : 
ibiock, "b" input or f 
output) 

CAP(a,b) ae 1p | farads ;Capacitors ("a"' gives 1 
| ‘ block, 7b" input or ' 

{ 

' | output) 

| 
| XL(a,b) C | P ohms Inductive reactance | 
i (same indexing as above) | 

XC(a,b) C i P | ohms ;Capacitive reactance | 
i(same indexing as above) + 

ZONE C ohms ‘Tee Leg one impedance 

| | 


Leg two impedance ' 


CO NOG Oe ERNE 


| 
} 
ZTHREE C ; ohms Leg: three impedance : 
5 { 

ZANT (L) C i ohms /Aerial input impedance | 
| ( ‘in frequency deviation | 
| | ! | | 

‘run 
{ } 
| 2HELIX(L) C | P ohms \Helix input impedance in } 
i ‘frequency deviation run | 
| ; q a 
| ZRES(L) | C ! P ' ohms ‘Bandwidth circuit [ 
| | | ,impedance in frequency | 
i deviation run 
| | 
| ZTEECL) eS ue iP | ohms \Tee input impedance in 
(frequency deviation run | 
ZTOWR(L) | C P | ohms jAerial and Coth@ae2) 
,impedance in frequency 
' ‘deviation | 
i i 
| onow | R ; P als Q of aerial circuit 
| | 
i { 


fe 
D 
ip) 
= 
Cc 
jan 
P- 
2} 
oo 
= 
le) 
n 
n 
5 
i) 
wn 
= 
1) 
ct 
ish) 
2 
(Pe) 
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Aerial self-resonant 
frequency 


SELRES 


XFMR1(L) 


} Transformer one unit, 
primary impedance, 
frequency run 


XFMR2(L) c ohms Transformer two units, 
} primary impedance, 


frequency run 


} 

| ohms Transformer three units, | 
primary impedance, 
| 

i 

| 


frequency run 


XFMR4(L) C ohms Transformer four units, 
primary impedance, 


XFMR3(L) C 
frequency run 


CURANT (L) C amperes; Aerial Current during 
(enol response 


run 


Current through ZTOWER(L) 
frequency run 


CURTWR(L) C amperes 


CURHEL (L) Helix current, frequency 


run 


amperes 


Bandwidth resistor 
current, frequency run 


CURRES (L) 


amperes 


Tee network current, 
frequency run 


CURTEE(L) amperes 


Transformer one unit, 
primary current, 
frequency run 


CURXL (L) 


CURX2(L) C 


a ad two units, 
primary current, 
frequency run 


| CURX3(L) € amperes} Transformer three units, 
primary current, 
frequency run 
CURX4(L) C amperes; Transformer four units, 


| 
t 
| 
| 
primary current, 
frequency run 


Aerial, forward, 
ABCD matrix 


Akon.) 


Helix, focward, 
ABCD matrix 


B(c,d) 


Bandwidth resistor, 
forward, ABCD matrix 


Ctexd) 


Ht) 


FS ONS A, CT TRIG ENO ET A See" SN “Nr 


ens nt AEN PN RR, 


WCced) 
Et egd) 


ECe,d) 


G(c,d) 
O(c,d) 


P(c,d) 


R(c,d) 


Ge. cd) 


U(c,d) 


V(c,d) 


W(c,d) 


X(c,d) 


EP(c,d) 


RP (c,d) 


UP(c,d) 


WP(c,d) 


CURRAD(L) 


PWROUT(L) 


ZTANK (a) 


Q 


QO 


Q 


rn PERE Semen AN ES SA SP ENR A RCE AERTS rt NY AAS te TET OO O° 


a a ree FS 
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a ee er eee 


Tee, forward, ABCD matri 


Transformer one, forward, 
ABCD matrix | 


- Transformer one, reverse 
ABCD matrix 


- Tee reverse, ABCD matrix 


Helix reverse, ABCD 
matrix 


Aerial reverse, ABCD 
matrix 


- Transformer two, forward, 
ABCD matrix 


- Transformer two, reverse, 
ABCD matrix 


- Transformer three, 
forward, ABCD matrix 


~ Transformer three, 
| reverse, ABCD matrix 
{ 


~ . Transformer four, 
| forward, ABCD matrix 


- | Transformer four, 
‘reverse, ABCD matrix 


- ‘Amplifier one, forward, 
| ABCD matrix 


- ‘Amplifier two, forward, 
_ABCD matrix 


- “Amplifier three, 
| forward, ABCD matrix 


= | Amplifier four, forward, 
| ABCD matrix 


amperes! Current through radia- 

Eee resistance, 

, frequency run 

watts ieenen across radiation 
| resistance, freq. tum 


ohms eerie secondary 
impedance with power 
amplifiers connected 
(a goes 1-4) 


V. Discussion 

The construction of an antenna system for a high power Very Low 
Frequency transmitter station is an extensive undertaking. The 
physical structure is usually very large (but electrically small), and 
is a compromise from the very beginning. .Aerial deficiencies cannot be 
completely overcome, but by a prudent choice of coupling and matching 
system full use may be made of the finally accepted system design. 

In the past much use has been made of basic theoretical 
studies followed by modeling studies of the most promising designs. 
Useful information from these tests has been accurate enough to be 
used in the final design and also to develop the necessary coupling 
and matching networks. 

A network analysis for the complete system is made first at 
the limits of. the operating frequency range to determine extremes for 
tunable components. A more detailed analysis is then carried out for 
the operating frequency. If the transmitted waveform involves fre- 
quency deviation (such as with FSK), the antenna system is analyzed 
at the frequencies corresponding to the maximum deviation encountered. 
Within this range intermediate values may be obtained by reference to 
a set of universal curves for the type of circuit involved. 

If the antenna system at Northwest Cape is used as an example, 
the circuit may be analyzed at resonance in terms of input and output 
parameters at the various network blocks in a relatively short time. 

A more detailed analysis, such as computation of voltage, current, and 
power dissipated for each of the individual elements, will take much 
longer. The complexity increases if the frequency is allowed to vary 


from resonance, such that a complete analysis becomes very time 
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consuming. The choice Of frequency steps when operating off 
resonance is important: is five cycles satisfactory or should 
increments of one-half cycle be used? 

Parameter values in the field will vary from those measured 
in model studies or found from a purely theoretical analysis. If 
the antenna system is designed with such changes in mind (provisions 
for shunt capacitance, additional load inductance, variation in loss 
resistance, etc.), then the final adjustement in the field will be 
more easily. accomplished. To.be prepared for such variations a good 
theoretical analysis is necessary, which requires a rather large 
amount of time. If a major parameter such as the radiation resistance 
is found to have a different value in the field, however, the design 
analysis must be modified in whole, since it is no longer correct, 

It is at this point that a computer program for network 
analysis can be extremely beneficial. A number of such programs have 
been developed in the last several years by IBM and others. The 
programs vary in their usefulness, since they are usually very 
general in design. Within the last year. the Jet. Propulsion Laboratory 
has introduced two programs for antenna. design, but these are applicable 
only for specialized antennas at much higher frequencies. While most 
of these programs could be of some use, more beneficial results could 
be obtained from a program tailored to the VLF antenna system in 
particular. The author is aware of only one program other than VLFANT 
which is so designed, and this program (part of a larger study) ap- 
parently has never been publicly dve@liosedivae " 

VLFANT is a general program, tailored to meet the special net- 


works encountered in a VLF antenna system. Once the initial parameters 
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have been defined (such as top-hat capacitance, radiated power, 
operating frequency, etc.), settings for tunable components may be 
determined. The frequency may then be varied at will, using 
increments of any desired size. Results are printed out and plotted 
to give added insight into system response if desired. Tuning ranges 
may be determined for variable components by choosing operating 
frequencies at the ends of the VLF band. Impedances, voltages, 
currents, phase angles, and so forth may be easily found. Any 
parameter may be modified and a new set of results will be available 
with only a few minutes of computer time. 

One area not covered in VLFANT is the operation of the Class C 
power amplifiers normally encountered in high power transmitter. In 
the analysis the power amplifiers are considered to be constant 
voltage generators. This is a reasonable assumption at sade resonant 
frequency where the antenna system load is designed to be a pure 
resistance, but off resonance the load will vary and include a 
reactance term, and the constant voltage assumption will no longer 
be valid. A thorough analysis is difficult, since the effective 
plate resistance (treated as the generator internal resistance) is 
a function of various parameters including the load impedance magnitude 
and phase angle. 

One empirical study of this problem has been made, where an 
actual circuit was constructed and various measurements made and 
plotted. While this gives a good indication of Class C operation in 
general, it is good only for the tube and circuit tested, and does 


not yield a theoretical generalization of the eamsienee. 
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One student at the Naval Postgraduate School has written a 
computer program involving the use of Class C amplifiers as part of 
his thesis work on network synthesis using active a eae Basic 
tube parameters (plate resistance, output capacitance, amplification 
factor, etc.) and maximum ratings (plate dissipation, grid drive, 
plate voltage, bias levels, etc.) are first entered for the tube 
under study. Then the minimum output power, tube efficiency and 
bandwidth desired at the operating frequency are entered. The 
program computes the necessary operating parameters to meet the 
required conditions, and indicates which requirement is critical if 
the desired conditions cannot be met. An LC tank circuit is the 
tube load and a value for the necessary inductance is given for 
desired operation. 

While this program is not immediately usable with VLFANT, 
it appears to be an excellent starting point for the inclusion of 
Class C amplifier operation. A combination of the Class C amplifier- 
VLFANT programs should give all information necessary for a complete 
VLF transmitter-antenna system analysis. 

A second area of interest is the transient response of the 
antenna system when fed with an FSK input. The ABCD matrices developed 
in VLFANT may be used to form a transfer function in a Laplace analysis 
of the circuit response, but due to limitations in time such a study 
was not undertaken in this thesis. It should be noted that in the 
circuit analysis reported in DECO Interim Report No. 34-C-6 an FSK 
signal using a modulation index of one produced less transient response 
degradation than did a similar signal with an index of one-half. In 


no case, however, did the transient distortion go to zero, so strict 
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attention should be paid to this subject since in systems operating 
near breakdown levels transients may cause severe damage to the system. 
With respect to the Northwest Cape VLF station used as the 
model in VLFANT, the output data was found to be in excellent agree- 
ment with that measured by the facility itself. The numerical 
examples used in the Basic Principles section of this thesis was 
taken from original design information supplied while the program 
was being written. After the transmitter installation was completed, 
thorough site measurements were made and some of the parameters were 
found to be noticeably different from what had been predicted. In 
particular, the loss resistance was lower than first thought, and a 
40 foot length of bus bar added some series inductance to the aerial. 
These variations were easily taken care of in VLFANT: the new data 
was entered into the program and with only a few minutes of computer 


time the new output data was available. 
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Wile CONCLUSIONS 

The design and construction of VLF transmitting antenna systems 
is an involved and expensive project. With the increased importance of 
VLF communications in the past few years efficient and reliable opera- 
tion of all components of such a system becomes even more important. 

The best analysis possible of a proposed system must be achieved, 
and all data relevant to existing system should be up-to-date, so that 
maximum overall efficiency is attained. 

Computer programs designed for use on high-speed digital com- 
puters help to achieve an accurate analysis of systems in a short time. 
VLFANT, designed expressly for VLF antenna systems, will provide infor- 
mation of relevance to this specialized field, and may be used from the 
beginning design stages through to a fully operational status. When 
tested with the Northwest Cape transmitter installation as a model, 
output data was in excellent agreement with data measurements made in 
the field. A large quantity of additional information, either too 
difficult or time consuming to actually measure, was also provided, 
which gives a deeper insight into the operation of the station. 

The digital computer program developed here has only scratched 
the surface of the many advantages available in using such a tool for 
the optimum design of VLF antenna systems. Some areas requiring ad- 
ditional study have been pointed out (Class C amplifier operation, 
transient response, etc.), and by an extended use of a program such 
as VLFANT, other such areas may be determined. 

Full application of such a program should be made, especially 
in the all important design stage, for as has been pointed out before 


there are two ways to design an antenna system: "(1) design the unit 
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properly, in a comfortable office during the day, or (2) slowly 


freeze, night after night, while trying to make a poor design orketae 
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APPENDIX A 
BANDWIDTH RESISTOR NETWORK: 
The "Q" of an aerial is determined by the largest reactance 
together with the total resistance present in the circuit, and for 


most VLF aerials may be expressed by: 


[Fe 
O28 en ane -+ C1) 
where eT = Quality factor for aerial 
x = Value of aerial static capacitance 
R = Resistance 


When speaking of the aerial alone, R is the radiation resistance; for 
the antenna system as a whole R is the sum of the radiation and loss 


resistances. See Figure 44. The corresponding bandwidth is given by: 


BW = “o (2) 
0 oie 
where BW = 3 db bandwidth in Hz 
fo = operating frequency in Hz 


Similar definitions for Q and bandwidth are applicable for 
series circuits at their resonant frequency. 

Networks in the complete antenna system are designed to be 
lossless, that is, they have extremely low loss resistances and have 
reactances much greater than the resistance present. When these 
networks are used for impedance matching at the operating frequency 
their reactance terms cancel out such that. no new reactance is added 
to the overall system, and being lossless no additional resistance 


is added. 


Hes 


R .dtecion 


Figure 44. Basic Aerial Circuit 
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Bandwidth 
Resistor 
Network 


Figure 45. 


| oh» _h_| 
3, 


Helix Network eee Network 


Basic Bandwidth Widening Circuit 


Figure 45 shows a helix impedance matching circuit con- 
nected to the basic aerial circuit. To the right of cut A the 
input impedance will appear as a pure resistance at the resonant 
frequency. There might be a tendency to apply Equations (1) and 
(2) at this input, but they no longer apply since the circuit is 
no longer a simple series circuit. To the left of cut A a new 
Bandwidth Resistor network has been connected. By varying 
inductance Ly the effect of resistance R, is changed. 

To properly evaluate the actual circuit Q when the Bandwidth 
Resistor is added to the network, a definite procedure must be used. 
In Figure 46 all series reactances in the helix-aerial circuit have 
been combined; from the design of the helix network. the resulting 
reactance will be capacitive. The series Bandwidth Resistor circuit 
has been put into its equivalent parallel form. 

Figure 47 shows the aerial RC' series network replaced by its 
equivalent parallel form - Roe and the parallel resistance of the 


bandwidth resistor has been moved to the far right hand end of the 


network. The resistances S - R,' may now be combined, and then the 


2 
resistance sum and ms combination put into series form, Rew’ * See 
Figure 48. 


The capacitive reactance undergoes little change during the 
above transformation, so C' may be put back into a form similar to 
the original circuit, as shown in Figure 49. Now a new Q for the 


circuit may be defined. 


Q = R wea) 


new 


The procedure described above is useful where an assumed 


value for the Bandwidth Resistor network coupling inductor is used. 


de 


Figure 46. Reactances Summed 


Figure 47. Parallel Representation 
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Figure 48. Series Form of Aerial Circuit 
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Figure 49. Resultant Widened Circuit 
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As L, is varied the resulting parallel resistance R,' will vary, 


4 2 


resulting in a range of values for OA 

This type of transformation above will hold for only one 
frequency since the reactance values effect the final value of Rew’ 
Additionally, the helix circuit was designed to match the original 
aerial resistance R, and the. desired input impedance will no longer 
be seen with ee in the circuit. If the circuit impedance is 
computed at other than the operating frequency, an even larger 
discrepancy will be noted. 

The question then is how to properly evaluate the effect of 
the Bandwidth Resistor network. This can be done by observing 
the current flowing into the series aerial circuit and through the 
radiation resistance. An ammeter is placed in series with the con- 
nection between the helix and aerial circuit, at the junction of 
Ly and Lo. Then when voltage is applied at the input to the overall 
circuit, Figure 45, current will divide into three paths: through 
Los L3 and Ly. 

In actual operation the voltage at, the. circuit input will be 
different with the Bandwidth Resistor..in the system since the overall 
impedance will draw more current through previous stages (such as a 


Tee impedance matching network) and the new system response may be 


properly evaluated. 
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SUBROUTINE ANITAMP { K, MONE, MASK,—HK ’ 


THIS SURPOUTINE WILL DETERMINE THE VARTATION 
IN CURRENT WITH FREQUENCY, USING THE VOLTAGE 
MEASIIREN AT EACH OF THE NONES WHERE CURRENT 
1S COMPUTEN, 
VOLTAGE AND CURRENT ARE MEASIIREM PLATE-TO-PLATF 
ANT WILL SF NETERAENEN BREFTRE THIS SURROTITINEG 
IS CALLEN. MATRIX SUBROUTINE WILL BE USED TO 
NETF2MENE VOLTAGES AND SUBROUTINE DFLTA THE 
EYPFOANCFES. GROUP ONE TNCLINES TRANSFORMER, 
TFE, RESESTAR. AND HEC TX CURRENTS. GeRonue Twn 
INCLIINES TOWER ( "ANTENNA" 9 , AERTALs 
ANT RANTATION RESTSTANCES CURRENTS AND POWFR 
OUTRUT IN WATTS. 
MOQF = ] USE 1 TRANSFORMER 

MADE = 2? H11SF 2 TRANSFORMERS 

MONF = 3 JSF 3 TRANSFORMERS 
MIOF = 4 USE 4&4 TRANSFENRMERS 
HK = % NEGLECT BANDWIOTH RESTSTOR 
Hk = } HSE RANNWIOTH RESTSTAR 

ASK =a} NO PRINTOUT 

MASK = 2 RRINT GROUP ONE 

MASK = 3 RRINT GROUP TWO 

MASK = G4 ORINT GROWPS NNFE fF TWO 
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ee FOUP AV ee REUA AV, UPT2.7Ade WHL2.7AV, WWI7.2) 
T7241), CIRTWRI241), CURHFLI2419, CURRES( 2411, 
CHRXL (241), CURX20 7414, CHIRXZ(2419, CURK4E241) 

Nf 761), PWROUITO2411 

(94)) 
LES. 4). CAP(5e34- RAS, RACT RAT, RAW, PT ia 
LTS. AMPS, XCe XbLe XKTOTALs VYTOTALs ZTOTAL 
QO741¥—6 NIVI261), TANT» 7HELIXs. 7RESs ZTEE. 7TOWR 
Rl» YEMR2, XFMR3, XEMRG 

ANT. CURTWR, CURHEL, CURRES., CURTFE 

Re Co Ds Fe Fe Ge Oe Py Re Te Ue Ve We Xo title VV 
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QP, IP, WP 
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HOMIE = rik 

FOFACY = FREO() 

CAME IMATRIY € Saracy, 9 

(alr wal ¢ Wo Bo ato il Yo Milalinte 

CueKi(ti = “P(1.1' 7 XEMRI(L) 
CHew(Ly = CURYI]I) 

FeR(?.19 = — FRTF.1) 

CALL MULT € FF, FO, Wwe 2e 25 Pe 1 ) 
VALTS#H2<et} = WH{4.1) === = me = oe 
BO Tor 21 

CuRK?(L)I = ie 7 X¥FMRO(LY) 
CUax(L) = Clb AFI 

DIP real ae = Ot >a) 

Caw MULT ( ie RE, Wil, Pe 2° Pe il 1 
VatTSE2.1) = Wwll.lt 

6A 9 Al 

(GUUS Seay WS (IP(J.1) 4 Y¥EMRAIL)D 
Cue = Gn vai) 

PN ZG ee ee NCR 

CALL MINT € Ve HIP, Wee Pe Fe De 1 ’ 
VOUTSE2.1) = wwll.d} 

BRA TN Pl 

CtuRxXalt) = WP{]T.1) 7 YEMRG(L) 
CUR¥YFt dy = CHRY4CI) 

AP(7.1) = = wR(?.1) 

CALL MULT € X¥. WE, WW, 276 Ae Pe 1 ) 
VO_TS¢@2<])) = wWwfd..) 

FIN YVEEIL YD = Mal ySteolth 2 vee t(li 
Fal d MIT € Re Whee ING De Pe Pe 1 ) 
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THTS SUBROUTINE WILL PLOT UP TO 9 NEPENNENT ARRAYS, Y1 THRU Y9, 
AGAINST THE [SNE PENDENT VARIABLE ®X®, ON THE PAGE PRINTER, THE 
AF POINTS EN THE X ARRAY WHEE E-PETERMENE- ¥ e 
ROTH THE X¥-Ax]S ANN Y-AXIS WILL BE PRINTED IF THEY ARE IN THE R 
THF GRAPH, THE MATIN PROGRAM MUST BE NIMENSTONED WITH THE NUMMY 
VARTASLFS ¥1L(1) THR ¥9(01) AND TITLELCL) THRU TITLE5(1). THESE 
VARTARLFS WELL RF USED WHEN LESS THAN 9 E€URVES ARE TO BF PRINTE 
SFE FXAMPLE AFI OW, THE X—-ARRAY MAY RFE SORTED INTO ASCENDING OR 
IF NECESSARY OR NFSTREN BY THE USER BY CALLING KK. THE PRONUCT 
OF THF NIIMRFR OF CURVES PLOTTEMD (PLUS 1) AND THF NUMBER NF POIN 
— AGH GHIRVE MHST-BE+FSS—THAN-BR-EQUAL— s——FAR- F Tt 
CURVES OF 25% POINTS FACHe. WILL HAVE. (341) * (250) = 1000. FN 
9 CURVES A MAXTMIIM NF 190 POINTS MAY RF USEN IN FACH. 
TITLING FOR THE CURVES WILL BE DONE THROUGH THE USE OF TITLFl T 
TITLES FACH CURVE NAME MAY BE 1P TO 28 CHARACTERS LANG. TIT 
USE THE FIOST >A SPACES FOR THE NAME OF CURVE #1, THe NEXT 12, 
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